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Physical interactions in quantum many-body systems are typically local: Individual constituents interact mainly 
with their few nearest neighbors. This locality of interactions is inherited by a decay of correlation functions, 
but also reflected by scaling laws of a quite profound quantity: The entanglement entropy of ground states. This 
entropy of the reduced state of a subregion often merely grows like the boundary area of the subregion, and not like 
its volume, in sharp contrast with an expected extensive behavior. Such "area laws" for the entanglement entropy 
and related quantities have received considerable attention in recent years. They emerge in several seemingly 
unrelated fields, in the context of black hole physics, quantum information science, and quantum many-body 
physics where they have important implications on the numerical simulation of lattice models. 
In this Colloquium we review the current status of area laws in these fields. Center stage is taken by rigorous 
results on lattice models in one and higher spatial dimensions. The differences and similarities between bosonic 
and fermionic models are stressed, area laws are related to the velocity of information propagation in quantum 
lattice models, and disordered systems, non-equilibrium situations, and topological entanglement entropies are 
discussed. These questions are considered in classical and quantum systems, in their ground and thermal states, 
for a variety of coiTelation measures. A significant proportion of the article is devoted to the clear and quantitative 
connection between the entanglement content of states and the possibility of their efficient numerical simulation. 
We discuss matrix-product states, higher-dimensional analogues, and variational sets from entanglement renor- 
malization and conclude by highlighting the implications of area laws on quantifying the effective degrees of 
freedom that need to be considered in simulations of quantum states. 
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I. INTRODUCTION 

In classical physics concepts of entropy quantify the extent 
to which we are uncertain about the exact state of a physical 
system at hand or, in other words, the amount of information 
that is lacking to identify the microstate of a system from all 
possibilities compatible with the macrostate of the system. If 
we are not quite sure what microstate of a system to expect, 
notions of entropy will reflect this lack of knowledge. Ran- 
domness, after all, is always and necessarily related to igno- 
rance about the state. 

In quantum mechanics positive entropies may arise even 
without an objective lack of information. To see this, let us 
consider a quantum lattice systems (see e.g.. Fig. 1) as an ex- 
ample for a quantum many-body system where each of the 
vertices i of the lattice L is associated with an individual quan- 
tum system. This quantum many-body system is thought to 
be in its non-degenerate pure ground state p = {tp\ at zero 
temperature which has vanishing von-Neumann entropy 

S{p) = -tr[plog2p]. 

Let us now distinguish a region of this quantum lattice system, 
denoting its sites with the set / and all other sites with O ~ 
L\I. If we consider the reduced state pi — Xxq [p] of the sites 
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FIG. 1 A lattice L with a distinguished set / C -L (shaded area). 
Vertices depict the boundary dl of / with surface area s(J) = \dl\. 



of the region /, the state will not be pure in general and will 
have a non-vanishing von-Neumann entropy S{pi). ' 

In contrast to thermal states this entropy does not originate 
from a lack of knowledge about the microstate of the sys- 
tem. Even at zero temperature we will encounter a non-zero 
entropy! This entropy arises because of a very fundamental 
property of quantum mechanics: Entanglement. This quite in- 
triguing trait of quantum mechanics gives rise to correlations 
even in situations where the randomness cannot be traced back 
to a mere lack of knowledge. The mentioned quantity, the en- 
tropy of a subregion is called entanglement entropy or geomet- 
ric entropy and, in quantum information, entropy of entangle- 
ment, which represents an operationally defined entanglement 
measure for pure states (for recent reviews see refs! '^^ ! '^ ^. 

In the context of quantum field theory, questions of scal- 
ing of entanglement entropies in the size of / have some tra- 
dition. Seminal work on the geometric entropy of the free 
Klein-Gordon field^^ and subsequent work on conformal 
field theories^^ '^^^^^ was driven in part by the intriguing 
sugges ted conn ection to the Bekenstein-Hawking black hole 
entropjEElIini. 

In recent years, studies of properties of the entanglement 
entr opy in thi s sense have enjoyed a revival initiated in 
j.gfg j7|i7i | i72 | 223[ Importantly, this renewed activity is benefit- 
ting from the new perspectives and ideas of quantum informa- 
tion theory, and from the realisation of their significance for 
the understanding of numerical methods and especially their 
efficiency for desciibing quantum many-body physics. Quan- 
tum information theory also provides novel conceptual and 
mathematical techniques for determining properties of the ge- 
ometiic entropy analytically. 

At the heart of these studies are questions like: What role do 
genuine quantum correlations — entanglement — play in quan- 
tum many-body systems? Typically, in such investigations, 
one abstracts to a large extent from the microscopic specifics 
of the system: Quite in the spirit of studies of critical phe- 



Of interest are also other entropies, such as the Renyi entropies, Sa{p) = 
(1 — a)~^ log2 tr[p"] with a > 0. For q \, 1 the usual von-Neumann 
entropy is recovered. In particular in the context of simulatability, Renyi 
entropies for arbitrary a play an important role. 



nomena, one thinks less of very detailed properties, but is 
rather interested in the scaling of the entanglement entropy 
when the distinguished region grows in size. In fact, for quan- 
tum chains, this scaling of entanglement as genuine quantum 
correlations — a prioii very different from the scaling of two- 
point correlation functions — ^reflects to a large extent the crit- 
ical behavior of the quantum many-body system, and shares 
some relationship to conformal charges. 

At first sight one might be tempted to think that the entropy 
of a distinguished region /, will always possess an extensive 
character Such a behavior is referred to as a volume scaling 
and is observed for thermal states. Intriguingly, for typical 
ground states, however, this is not at all what one encounters: 
Instead, one typically finds an area law, or an area law with 
a small (often logarithmic) correction: This means that if one 
distinguishes a region, the scaling of the entropy is merely 
linear in the boundary area of the region. The entanglement 
entropy is then said to fulfill an area law. It is the purpose of 
this article to review studies on area laws and the scaling of 
the entanglement entropy in a non- technical manner 

The main four motivations to approach this question 
(known to the authors) are as follows: 

• The holographic principle and black hole entropy: 

The historical motivation to study the entanglement or 
geometric entropy stems from considerations of black 
hole p hysics : It has been suggested in the seminal work 
of refs.l^EoIl that the area law of the geometric entropy 
for a discrete version of a massless free scalar field — 
then numerically found for an imaginary sphere in a ra- 
dial symmetry — could be related to the physics of black 
holes,^^^ in particular the Bekenstein-Hawking entropy 
of a black hole which is proportional to its bound- 
ary suiface. It has been muted that the holographic 
principle- — the conjecture that the information con- 
tained in a volume of space can be represented by a the- 
ory which lives in the boundary of that region — could 
be related to the area law behavior of the entanglement 
entropy in microscopic theories. 

• Distribution of quantum correlations in quantum 
many-body systems: Area laws also say something 
quite profound on how quantum correlations are 
distributed in ground states of local quantum many- 
body systems. Interactions in quantum many-body 
systems are typically local, which means that sys- 
tems interact only over a short distance with a finite 
number of neighbors. The emergence of an area 
law then provides support for the intuition that short 
ranged interactions require that quantum correlations 
between a distinguished region and its exterior are 
established via its boundary surface. That a strict 
area law emerges is by no means obvious from 
the decay of two-point correlators, as we will see. 
Quantum phase transitions are governed by quantum 
fluctuations at zero temperature, so it is more than 
plausible to observe signatures of criticality on the 
level of entanglement and quantum correlations. This 
situation is now particularly clear in one-dimensional 
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systems I''I7I13I36I38I44I46I69I72I73I83I85I92I111I132I134I138I146I147I22: 

but pro gress has also been made in hi gher-dimensional 
sy3tcmsP^^'^'^^''^°''^°''^°''^^''^''^''^'^3^' with rigorous 
area laws specifically for quasi-free bosonicP^I^IliMi and 
fermionic^* systems, as well as in disordered 
system^Sni 

• Complexity of quantum many-body systems and 
their simulation: One of the key motivations for study- 
ing area laws stems from a quite practical context: The 
numerical simulation of quantum many-body systems. 
In fact, if there is little entanglement in a ground state 
of a many-body systems, one might suspect on intuitive 
grounds that one can describe this ground state with 
relatively few parameters. More specifically, for one- 
dimensional systems, one would expect numerical al- 
gorithms like the powerf ul density-matrix renormaliza- 
tion group method^^^(DMR.G) to perform well if the 
ground state contains a small amount of entanglement. 
This suspicion can in fact be made rigorous'" '^^ ^'^ 
as it turns out that the scaling of entanglement speci- 
fies how well a given state can be approximated by a 
matrix-product state**'* ''-"^ as generated in DMRG. It is 
hence not the decay behavior of correlation functions as 
such that matters here, but in fact the scaling of entan- 
glement. 

• Topological entanglement entropy: The topological 
entanglement entropy is an indicator of topological 
order a new kind of order in quantum many- 
body systems that cannot be described by local or- 
der parameterPsEomillllml Lattice models having 
a non-vanishing topological entanglement entropy may 
be seen as lattice instances of topological quantum field 
theories. Here a global feature is detected by means of 
the scaling of geometric entropies. 

In this Colloquium we do not have sufficient space to give 
an account of all known derivations of area laws for the entan- 
glement entropy. However, we will try not to merely remain 
at a superficial level and only state results, but will explain 
a number of key techniques and arguments. When we label 
main statements as "theorems" this is done to highlight their 
special role, to make it easier to follow the line of reason- 
ing. For details of arguments and proofs, often technically 
involved, we refer the reader to the original work. The reason 
for the technicality of proofs originates from the type of ques- 
tion that is posed: To distinguish a region of a lattice breaks 
the translational symmetry of the problem - even in a transla- 
tionally invariant setting. While numerical studies are some- 
times easier to come by, analytical argument can be techni- 
cally involved, even for quasi-free models. In this article, we 
discuss the study of entanglement entropy primarily (i) from 
the viewpoint of quantum information theory, (ii) with an em- 
phasis on rigorous and analytical results, and (iii) the implica- 
tions on the efficiency of numerical simulation. 



<ll. LOCAL HAMILTONIANS AND AREA LAWS 

Throughout this article, we will consider quantum many- 
body systems on a lattice. Such quantum lattice systems are 
ubiquitous in the condensed matter context^^^ where they play 
a key role in obtaining an understanding of material proper- 
ties from a microscopic basis. Lattices systems are also of 
considerable importance in the study of quantum field theo- 
ries where a lattice provides a natural ultra-violet cut-off and 
facilitates numerical simulations of quantum fields'^''. One 
could think, e.g., of systems of strongly correlated electron 
systems or lattice vibrations of a crystal lattice. With the ad- 
vent of research on cold atoms in optical lattices, quantum 
lattice systems can also be prepared in laboratory conditions 
with an unprecedented degree of control^2i 

We will consider — at least in parts of this article — general 
lattices. Each vertex of the lattice is associated with a quan- 
tum system, such as a spin, a bosonic or a fermionic sys- 
tem. It is convenient to think of this lattice as a simple graph 
G = (L, E) with vertices L, and the edge set E labeling 
neighborhood relations. G could be the graph representing 
a one-dimensional chain with periodic boundary conditions, 
and in fact a good proportion of this article will deal with 
such quantum chains. For later purposes, it will be conve- 
nient to think in terms of such a slightly more general picture, 
however The Hilbert space of the total many-body system is 
then the tensor product 

where Hj is the Hilbert space associated with the physical 
system on lattice site j. On such a lattice, one has dist(j, k) 
for j,k G L as the natural graph theoretical distance which is 
the length of the shortest path connecting j and k. For a cubic 
lattice of dimension T) with periodic boundary conditions, in 
turn, dist(j, k) — X^dLi \jd — kd], where the components of 
j,k E L are taken modulo the base length of the cubic lattice. 

We will be concerned largely with local Hamiltonians on 
lattices. This means that the physical system associated with 
a specific lattice site will interact only with its neighbors and 
not with all sites of the lattice. The total Hamiltonian can 
hence be written as 

XCL 

where Hx has a compact support X, independent of the sys- 
tem size, that is the number of lattice sites denoted by \L\. 

The boundary surface area s{I) of a distinguished region 
/ of the lattice L can be defined in a very natural fashion on 
such a graph as the cardinality of the set of boundary points 

dl ^ {i e I : there is a j e L\I with dist{i,j) = 1} , (1) 

so ,s(/) = \dl\, see Fig.[T| Throughout the article, unless de- 
fined specifically otherwise, we will say that the entanglement 
entropy satisfies an area law if 

S{pi) = 0(s(/)). 
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This means that the entropy of the reduced state pi scales at 
most as the boundary area of the region /. 

Before we dive into the details of known results on area 
laws in quantum many-body systems, let us appreciate how 
unusual it is for a quantum state to satisfy an area law. In 
fact, a quantum state picked at random will exhibit a very 
different scaling behavior. If one has a lattice system with 
d-dimensional constituents and divides it into a subsystem 
I C L and the complement O = L\I, then one may consider 
the expected entanglement entropy of / for the natural c hoice, 
the unitarily invariant Haar measure. One find^^^IIllol] 

dl^l-ioi 

E[5(,,)]>|/|log,(d)-^^. 

That is, asymptotically, the typical entropy of a subsystem 
is almost maximal, and hence linear in the number of con- 
stituents |/|. Hence a "typical" quantum state will asymptot- 
ically satisfy a volume law, and not an area law. As we will 
see that area laws are common for ground states of quantum 
many-body systems, we find that in this sense, ground states 
are very non-generic. This fact is heavily exploited in numer- 
ical approaches to study ground states of strongly correlated 
many-body systems: One does not have to vary over all quan- 
tum states in variational approaches, but merely over a much 
smaller set of states that are good candidates of approximat- 
ing ground states of local Hamiltonians well, that is states that 
satisfy an area law. 

III. ONE-DIMENSIONAL SYSTEMS 

Most known results on area laws refer to one-dimensional 
chains such as harmonic or spin chains. This emphasis 
is no surprise: After all, a number of physical ideas — like 
the Jordan-Wigner transformation — as well as mathematical 
methods — such as the theory of Toeplitz determinants and 
Fisher-Hartwig techniques — are specifically tailored to one- 
dimensional translationally invariant systems. 

If we distinguish a contiguous set of quantum systems of 
a chain, a block I = {1, . . . ,n} the boundary of the block 
consists of only one (two) site(s) for open (periodic) boundary 
conditions. An area law then clearly means that the entropy is 
upper bounded by a constant independent of the block size n 
and the lattice size \L\, i.e., 

S{pi) = 0(1). (2) 

We will see that in quantum chains, a very clear picture 
emerges concerning the scaling of the entanglement entropy. 
Whether an area law holds or not, will largely depend on 
whether the system is at a quantum critical point or not. We 
will summarize what is known in one-dimensional systems at 
the end of the detailed discussion of quantum chains, starting 
with bosonic harmonic chains. 

A. Bosonic harmonic chain 

Bosonic harmonic quantum systems, as well as fermionic 
models and their quantum spin chain counterparts like the XY 



model, play a seminal role in the study of quantum many -body 
systems. Harmonic lattice systems model discrete versions of 
Klein-Gordon fields, vibrational modes of crystal lattices or 
of trapped ions and serve generally as lowest order approxi- 
mations to anharmonic systems. The fact that they are inte- 
grable renders even sophisticated questions like the scaling of 
the geometric entropy in instances amenable to fully analyt- 
ical study, even in higher spatial dimensions. In fact, in the 
latter case these so-called quasi-free models are the only set- 
tings that allow for rigorous analytical results so far Hence, 
they do form the central object of consideration to explore 
what should be expected concerning general scaling laws. 
The Hamiltonian for a harmonic lattice L is given by 

-ff = ^ ^ [PiPijP] + XiXijXj), (3) 

where X,P £ rI^I^^I^I are real, symmetric and positive ma- 
trices determining the coupling structure of the systems. The 
canonical operators Xi,pi satisfy the canonical commutation 
relations [xj , pfe] = In terms of the bosonic annihilation 
operators bj = [xj + ipj)/\/2 the Hamiltonian eq. (jsj reads 

H=\Y. (^I^^J^J + + ^^B^^]bJ + b\B,,jb]), (4) 

where A ^ {X + P) /2, B = (X - P)/2. Ground and ther- 
mal states of the above Hamiltonian are fully characterized 
by the second moments of the canonical operators, while first 
moments vanish'^ (entanglement properties of the state are in- 
variant under changes of first moments anyway). The second 
moments define the covariance matrix 

Tjj- = Tj}) = {riTj) + (r^Ti), (5) 

where r — (xi, . . . ,x\i^\,pi, . . . is the vector of canon- 

ical operators. The covariance matrix of the ground state of 
eq. (|3]l is given by F = F^; ® Fp, where 

Fp=Xl/2(xl/2pxl/2)-l/2^1/2 

and Vx = Fp^, see refs.l^S^E] On the level of covari- 
ance matrices unitary operations express themselves as sym- 
plectic transformations S that preserve the commutation re- 
lations ak,i = i[rk,ri], i.e., ScrS'^ = a. Importantly, 
Williamson 's Theorem states that for any strictly positive ma- 
trix A e ^^N-x.2N j-jjgj-g exist a symplectic transformation S 
such that SAS'^ = D, where D is a diagonal matrix with the 
same spectrum as the positive square roots of (iaA)^. The 
eigenvalues di of D are called the symplectic eigenvalues of 
A. 

Now, what is the entanglement content of the ground state? 
To answer this we need to define entanglement measures and 
compute them in terms of the properties of the covariance ma- 
trix. The first of these is of course the entropy of entangle- 
ment. Williamson's theorem shows that any function of a state 
that is unitarily invariant is fully determined by the symplectic 
eigenvalues. Notably, the entropy of a Gaussian state p with 
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symplectic eigenvalues di, . . . ,djv of the co variance matrix 
of p is given by 

siP) - E log. ^ - V '^^^ V) ■ 

A key ingredient in the analytical work is another full en- 
tanglement measure that was defined in quantum information 
theory, the logarithmic negativity^^'''*^^^^^^^^'*^*''. It is de- 
fined as 

EN{p,I) = \og.,\\p^'\\,, 

where ||v4||i — tr[(A^^A)^/^] is the trace norm and p^' is the 
partial transpose of p with respect to the interior /. The 
partial transpose w.r.t. the second subsystem is defined as 
(|i)(fc| ® \ = \i){k\ (g) \l){j\ and extended by linear- 

ity. On the level of covariance matrices the partial transpose is 
partial time reversal, i.e. pi i— >■ —pi if i e / while Xi remains 
invariant. Then for p with covariance matrix F = Fp 
we find that p^' has covariance matrix F' = Fj; (FTpF), 
where the diagonal matrix F has entries Fij — ±i5ij, de- 
pending on whether a coordinate is in / or O: Then one finds 
for a state with covariance matrix F = T^^Tp the logarithmic 
negativity^— 

1 

En{p,I) = -Elog2niax{l,Afe(F;iFF-iF)}, 

k = l 

where the {Xk} denote the eigenvalues. The logarithmic neg- 
ativity has two key features. Mathematically, the importance 
of En{p, I) is due to 

ENip,I)>Sipi) (6) 

which holds for all pure states p. This upper bound for the 
entanglement entropy is simpler to compute as one does not 
have to look at spectra of reductions pi but of the full sys- 
tem. This renders a study of area laws possible even in higher 
dimensional systems. Secondly, in contrast to the entropy 
of entanglement, the negativity is also an entanglement mea- 
sure for mixed states, such as thermal states and provides 
an upper bo und on othe r important measures of mixed state 
entanglemenPEElEni. 

All of the above holds for general lattices L but for the 
moment we will focus on the one-dimensional setting, that 
is L = {1, . . . , N} where N is even to allow us to consider 
the symmetrically bisected chain I — {1, . . . , N/2\ with pe- 
riodic boundary conditions and P = 1. We concentrate on 
the ground state and discuss thermal states later.^ It is worth 
noting that in higher spatial dimension the natural analog of 
this setting, the half-space, is of some importance as it allows 
for a reduction of the problem in question to the 1-D case 



^ We do not discuss the entanglement properties in excited states liere as this 
area has not been explored in detail so fai^i9«22L. 



discussed her#^. Furthermore, the scaling behavior of the en- 
tanglement of the half-chain has direct consequences on the 
availability of efficient representations of the state by means 
of matrix-product states as will be discussed in some detail 
later on in this article. For a general nearest-neighbor cou- 
pling this means that X is the circulant matrix, 

X = circ(a,6,0,...,0,6), (7) 

as a consequence of translational invariance. b specifies the 
coupling strength, a defines the on-site term, \min{X) = a — 
2|6|, i.e., positivity demands a > 2\b\, and the energy gap 

1 /2 

above the ground state is given by AE — \J^^^ {XP) = {a — 
2|6|)^/^. For the logarithmic negativity of the symmetrically 
bisected half-chain we 

Theorem 1 (Exact negativity of the half-chain) Consider a 
Hamiltonian of a harmonic chain on L = {1, . . . , N} with 
periodic boundary conditions, P = 1, and nearest-neighbor 
interactions as in eq. Q. Then the entanglement entropy of 
the symmetrically bisected chain and the logarithmic negativ- 
ity satisfy 

(8) 

1/2 

where \\ ■ \\ is the operator norm and AE — XJ-^ {X). 

The quantity \\X\\ will later be related to the speed of sound 
in the system. This expression for the block entanglement 
quantified with respect to the negativity is exact and no ap- 
proximation. This was to the knowledge of the authors a first 
rigorous area law for a lattice system, complementing earlier 
seminal work for fields^''. Remarkably, this expression is en- 
tirely independent of N, the system size. The most important 
observation here is that an area law holds, which can be ex- 
pressed in terms of the spectral gap in the system: Whenever 
the system is non-critical in the sense that the energy gap AE 
satisfies AE > c > with a system size independent con- 
stant c, a one-dimensional area law will hold. The above link 
of entanglement entropy and spectral gap in the system can 
be established in much more generality and we will delay this 
discussion to a later subsection. 

The argument leading to Theorem[T]is involved, and for de- 
tails we refer to ref.^. The interesting aspect of this proof is 
that the spectrum of the half chain can not be obtained ana- 
lytically, thus not allowing for a direct computation of the en- 
tanglement content. Instead, it is the particular combination 
of spectral values of the partial transpose entering En{p,I) 
itself that can be explicitly computed. The proof makes heavy 
use of the symmetry of the problem, namely the invariance 
under a flip of the two half chains. 

This result suggested that the locality of the interaction in 
the gapped model is inherited by the locality of entanglement, 
a picture that was also later confirmed in more generality. 
Note that the above bound is a particularly tight one, and that 
it may well suggest what prefactor in terms of the energy gap 
and speed of sound one might expect in general area laws, as 
we will discuss later. 
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Let us now consider an important model for which the en- 
ergy gap vanishes in the thermodynamical Hmit N oo: 
Taking a = m? + 2N'^, h = —N"^, identifying lattice sites 
by j = xN, and the canonical operators by Xi = N^^^^(j){x), 
Pi = N~^/^7r{x), one obtains the Klein-Gordon field Hamil- 
tonian 



2 



dx \ Tr'\x) + ( ^(l^ix) ] + m^(t>\x) 



(9) 



in the field limit TV -> cjo. (For a detailed discussion of the 
continuum limit for the Klein-Gordon field, see also refP^.) 
From the expression ([8]l for the entanglement, we immediately 
obtain 



1 



En{p, I) = -\ogA I 



1 



log2 



2N 



(10) 

This is a striking difference to the area laws that we have ob- 
served earlier, now the entanglement does not saturate but di- 
verges with the length of the half-chain. The behavior ob- 
served here will be mirrored by a similar logarithmic diver- 
gence in critical quantum spin chains and fermionic systems. 
This will be discussed in the following section. 



are quadratic in fermionic operators f i and fj , 



(11) 

may be treated by similar analytical techniques and follow 
similar intuition to the bosonic case. In eq. (fTTb, to ensure 
Hermiticity of the Hamiltonian, = A and lP'= —B must 
hold for the matrices A and B defining the coupling. The role 
of the canonical coordinates is taken by the Majomna oper- 
ators Xj ^ (/] + fj)/V2 and = i(/] - fj)/V2, while 
the role of symplectic transformations is taken by orthogonal 
transformations. The energy gap above the ground state is 
given by the smallest non-zero singular value of A + B. 

Note that, in contrast to the bosonic case, the ground state 
is 2l^l"™''(^+^) =: g-fold degenerate. We define the ground 
state expectation (•) = tr[-Po]/<Z> where Pq projects onto the 
ground state sector Then, as in the bosonic case, the ground 
state is fully characterized by two-point correlations embodied 
in the covariance matrix with entries 

where now r — {xi, . . . . . . ,P\l\) collects Majorana 

operators. One then finds 



B. Fermionic cliain and the XY model 

Following the initial work on bosonic models of ref.' sim- 
ilar questions were explored in fermionic systems and the as- 
sociated spin models. The numerical studies in refs.EHSlpre- 
sented a significant first step in this direction . Their key obser- 
vation, later confirmed rigorouslyS^^EMIIl using techniques 
that we will sketch in this section, is that the scaling of the en- 
tanglement entropy as a function of the block size appears to 
be related to the system being quantum critical or not. Again, 
for a gapped system, away from a quantum critical point, the 
entanglement entropy would saturate, i.e., an area law holds. 
In turn in all cases when the system was critical, the numerical 
study indicated that the entanglement entropy grows beyond 
all bounds. More specifically, it grows logarithmically with 
the block size. This behavior is also consistent with t he be- 
havior of geometric entropies in conformal field theorjBHHl 
which applies to the critical points of the models discussed in 
refs.!^-* 147 223 -pjjg intriguing aspect here is that being crit- 
ical or not is not only reflected by the scaling of expectation 
values of two-point correlators, but in fact by the ground state 
entanglement, so genuine quantum correlations. 

This section defines the setting, introduces the basic con- 
cepts required and outlines the rigorous results in more de- 
tail. Fermionic quasi-free models, that is Hamiltonians that 



^ Compare also the divergence of the entanglement entropy in collectively 
interacting chain^iJj. 



yT Q V =\A + B\ + {A + B), (12) 



where •+ indicates the Moore-Penrose generalized inverse of 
a matrix'^"*, i.e., for a unique ground state one simply has V — 
\A + B\^^{A + B). The entropy of a contiguous block / of 
fermions in the ground state can be expressed in terms of the 
singular values ak of the principle submatrix Vi of V . One 
finds S{pi) = Y.k fi^k), where 



fix) = -^-7r-iog2 



l-x 



1 



l0g2 



1 



(13) 

All the above holds for general lattices but for the moment we 
will turn to a discussion of L = {1, . . . , TV}. 

We have started the discussion on the level of fermionic 
operators to highlight the similarity to the bosonic case. It is 
important to note however, that these fermionic models share 
a close relationship to natural spin models in the 1-D setting. 
This is revealed by the Jordan Wigner transformation which 
relates fermionic operators with spin operators according to 



= 1 - 2///., 



= n(l-2./lA)/., (14) 



k=l 



where af,af,af denote the Pauli operators associated with 



site i & L. The fermionic model eq. ( 1 1 1 is hence equivalent 
to a spin model with short or long-range interactions. 

The most important model of this kind is the XY model 
with a transverse magnetic field, with nearest-neighbor in- 



teraction, Ai 



A, A, 



-1/2 if dist(i,j) = 1, and 



7 



B, 



-Bj^i =7/2 for dist(i, j) = 1. This gives rise to 



H 



1+7 



T y 11 



ieL 



(15) 

where denotes summation over nearest neighbors, 7 is 
the anisotropy parameter, and A an external magnetic field."* 
Once again, translational invariance of the model means that 
the spectrum can be readily computed by means of a discrete 
Fourier transform. One obtains 



Ek ^ ((A - cos(27r/c/iV))2 + 7^ am'^{2TTk/N)) 



1/2 



,iV. 



This is a well-known integrable 



for k = 1, 
modePmi. 

In the plane defined by (7, A) several critical lines can be 
identified: Along the lines |A| ~ 1 and on the line segment 
7 = 0, |A| < 1, the system is critical, limAr_j.oo AE{N) — 
0. For all other points in the (7, A) plane the exists a c > 
independent of N such that AE > c. The class of models 
with 7 = 1 are called Ising model. The most important case 
subsequently is the isotropic case of the XY model, then often 
referred to as XX model or isotropic XY model. This is the 
case when 7 = 0. The XX model is critical whenever | A| < 1. 
The XX model is equivalent to the Bose-Hubbard model in the 
limit of hard-core bosons, so the Bose-Hubbard model with 
the additional constraint that each site can be occupied by at 
most a single boson. 

Let us assume that we have a non-degenerate ground state, 
such that the entropy of entanglement S{pi) really quantifies 
the entanglement content. For the translation-invariant system 
at hand, the entries Vi,j = Vi-j of V are given by 



V, 



fe = l 



9k 



A + B 



(16) 



The entanglement properties of the model are encoded in the 
numbers g^. For N = \L\ — > 00, we can write 



1 



27T 



"* Note that the boundary conditions give rise to a (sometimes overlooked) 
subtlety here. For open boundary conditions in the fermionic model, the 
Jordan Wigner transformation relates the above fermionic model to the spin 
model in eq. (B) with open boundary conditions. For periodic boundary 
conditions, the term /j^/i is replaced by the operator (Ylj^^iC^fj fj ^ 
l))/jv/l- Hence, strictly speaking, the peiiodic fermionic model does not 
truly correspond to the periodic XY model*' . Importantly, the degeneracy 
of ground states is affected by this. For a degenerate ground s tate the 
entanglement offormatioj^, the relative entropy of entanglement^, the 
distillable entanglemen^^or the logarithmic negativity of the ground state 
sector are the appropriate entanglement quantifiers and no longer the 
entropy of entanglement. Only in the case that for large subsystems n one 
can almost certainly locally distinguish the finitely many different ground 
states, the entropy of entanglement for each of the degenerate ground states 
still gives the coiTect value for the entanglement of a subsystem. 



for |/| < N/2, where g : [0, 27r) C is called the symbol 
of V. Note that the Fermi surface is defined by the discon- 
tinuities of the symbol. In order to evaluate the entropy of 
a reduction S{pi), we merely have to know the singular val- 
ues of n X n-submatrices Vi = Vn of V, see eq. (13 1. For 
isotropic models, i.e., for B = Q, V then being symmetric, 
the singular values are the absolute values of the eigenvalues. 
In other words, in order to understand the correlation and en- 
tanglement structure of sub-blocks of such systems, one has 
to understand properties of matrices the entries of which are 
of the form Ti j = Ti^j. Such matrices are called Toeplitz 
matrices. An n x n Toeplitz matrix is entirely defined by the 
2n — 1 numbers T;, ^ = 1 — n, . . . , n — 1. 

The spectral values Ai(r„), . . . , Xn{Tn) of T„ are just the 
zeros of the characteristic polynomial 

n 

det(T„-Al) = [|(Afe(r„)-A), 

fe=i 

so in order to grasp the asymptotic behavior of the spectrum of 
Tn, it is sufficient to know the asymptotic behavior of this de- 
terminant expression.^ The mathematical theory of determi- 
nants of such Toeplitz matrices is very much developed. The 
Fisher-Hartwig Theorem provides exactly the tools to study 
the asymptotic behavior of Toeplitz determinants in terms of 
the symbol. Crudely speaking, what matters are the zeros of 
g and the jumps: Once g is written in what can be called a 
normal form, one can "read off" the asymptotic behavior of 
the sequence of Toeplitz determinants defined by this symbol. 
Note that the matrices Vn — zl take here the role of T„. The 
exact formulation of the Fisher-Hartwig Theorem is presented 
in the appendix. 

This machinery was used in refs.l^SlIISIIIIlto evaluate the 
asymptotic behavior of the block entropy for the critical XX 
model and other isotropic models. In the first paper intro- 
ducing this idea'-'^, in fact, there is a single jump from 1 to 
— 1 in the symbol defining the Toeplitz matrices (and no ze- 
ros), which gives rise to the prefactor of 1/3 of log2(n) in 
the formula for the entanglement entropy in the XX model. 
This prefactor — which emerges here rather as a consequence 
of mathematical properties of the symbol — is related to the 
conformal charge of the underlying conformal field theory. 
In more general isotropic models, as has been pointed out 
in ref.'^'^, the number of jumps determines the prefactor in 
the entanglement scaling. Hence in such quasi-free isotropic 
fermionic models, the connection between criticality and a 
logarithmic divergence is very transparent and clear: If there 
is no Fermi surface at all, and hence no jump in the symbol. 



' Once this quantity is known, one can evaluate the entropy by means of 

S{pi) = J2 /(l^fcl) = / l°gdet(Vn - zl). (17) 

f—' 27ri dz 

fc— 1 

Here, the integration path in the complex plane has been chosen to contain 
all eigenvalues A^; (Vn). The function is a continuation of /: We require 
that limE:_»o /^(z) = including the parameter e such that is 

analytic within the contour of integration. 
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the system will be gapped and hence non-critical. Then, the 
entropy will saturate to a constant. 

In contrast, in case there is a Fermi surface, this will lead 
to jumps in the symbol, and the system is critical. In any 
such case one will find a logarithmically divergent entangle- 
ment entropy. The prefactor is determined by the number of 
jumps. So more physically speaking, what matters is the num- 
ber of boundary points of the Fermi surface in the interval 
[0,27r). So — if one can say so in a simple one-dimensional 
system — the "topology of the Fermi surface determines the 
pref actor". Th is aspect will later be discussed in more detail. 
Refs.EHMEllfind the following: 

Theorem 2 (Critical quasi-free fermionic chains) 

Consider a family of quasi-free isotropic fermionic Hamil- 
tonians with periodic boundary conditions as in eq. ([77J 
with B — 0. Then, the entanglement entropy of a block of 
/ = {1, . . . , n} continuous spins scales as 

S{pi)^^\og^{n) + 0{l). 

^ > /s a constant that can be related to the number of jumps 
in the symbol (defined above). This applies, e.g., to the scaling 
of the entanglement entropy in the XX spin model, for which 

5(p/) = ^log2(n) + 0(l). 

The constant ^ is not to be mistaken for the conformal 
charge which will be discussed later These arguments cor- 
respond to the isotropic model with B — Q, where the Fisher- 
Hartwig machinery can be conveniently applied. In contrast, 
the anisotropic case, albeit innocent-looking, is overburdened 
with technicalities. Then, in order to compute the singular 
values of submatrices of V as in eq. ( [T2] i, it is no longer suffi- 
cient to consider Toeplitz matrices, but block Toeplitz matrices 
where the entries are conceived as 2 x 2 matrices. This setting 
has been studied in detail in ref.E^ in case of a non-critical 
anisotropic system, finding again a saturation of the entangle- 
ment entropy and in ref where the prefactor of the area law 
for the entanglement entropy in the gapped XX model was 
computed rigorously. Ref.l^Sl discusses also other Renyi en- 
tropies in this model. 

Using an idea that originates from the concept of single- 
copy entanglement all these technicalities may be avoided and 
we can prove that the entanglement entropy diverges at least 
logarithmically in case of a critical (anisotropic) Ising model. 
The il notation just means that there is asymptotically a lower 
bound with this behavior.^ 

Theorem 3 (Divergence for the critical Ising model) The 

entanglement entropy in the critical Ising model scales as 



The starting point leading to this result from ref.E3is a lower 
bound in the operator norm of pj leading to 



Sipi) = Q{\og^n). 



(18) 



l0g2 



> 



log2 det((l 
;^log2|det(K)| 



Vn)/2) 



(19) 



This makes a big difference: We now no longer need the 
singular values of Vn (which would lead to an enormously 
complicated block Toeplitz expression, for a case for which 
the Fisher-Hartwig conjecture has not yet been proven). 
Instead — as the absolute value of the determinant is just as 
well the product of the absolute values of the eigenvalues as of 
the singular values — we can use the ordinary Fisher-Hartwig 
machinery to get an asymptotic handle on eigenvalues. For the 
critical Ising model, we can again find an explicit factorization 
of the Fisher-Hartwig-symbol, in terms of a function reflecting 
a single discontinuity and an analytical function. Using again 
a proven instance of the Fisher-Hartwig conjecture'^^ — albeit 
a different one than used in the case of an isotropic model — 
one finds the bound as in eq. ( 18 i; for details see ref.^^. The 
entanglement in two blocks of the critical Ising model has 
been studied in ref.l^. 

Another useful starting point to obtain bounds to entan- 
glement entropies in fermionic systems is to make use of 
quadratic bounds to the entropy function: Such quadratic 
bounds immediately translate to a bound to the function / in 
the expression of the entropy of a fermionic state in terms of 
the covariance matrix as in eq. ( [T3| ), as 



(20) 



This immediately translates to a bound of the form tr[(l — 
ViVfy^^] > S{pi) > tr[l - ViVf], where Vi is the sub- 
matrix of V associated with the interior 7 --. These bounds 
have also been exploited in the higher-dimensional analysis in 
ref^li. 

A method to obtain area laws in particular for symmetri- 
cally bisected quantum chains is the so-called method of cor- 
ner transfer matrices. This method has first been used in ref.l^ 
for the computation of the entanglement entropy, using ideas 
going back to ref.^^. The infinite sum of ref.-''' could be per- 
formed in ref.l^, giving also rise to a formula for the entan- 
glement entropy in the XX model. This idea has also applied 
to further models in ref!^-'". 

To conclude the discussion of critical quasi-free fermionic 
models let us note that the correspondence of being critical 
(gapped) and having a logarithmically divergent (saturating) 
entanglement entropy holds true for local systems only. If 
one allows for long-ranged interactions, then one can indeed 
find gapped, non-critical models that exhibit a logarithmically 
divergent entanglement entropy : 

Theorem 4 (Gapped model with long-range interactions) 

There exist models with long-range interactions, the coupling 
strength being bounded by r/dist(j,k) for some constant 
r > 0, such that for some constant ^ > 



«/(") = !^(s(n))if3C>0,no :Vn>no : \Cgin)\ < \f(n)\. 



^(p/)=Clog2(n)+0(l). 
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Hence, being gapped — albeit having power-law 
correlations — does not necessarily imply an area law. If 
one allows for long-range interactions (and a fractal structure 
of the Fermi-surface), one can show that one can even 
approa ch ar bitrarily well a volume law for the entanglement 
entropjEsEH Interestingly, states that are defined by quantum 
expanders can have exponentially decaying correlatio ns and 
still have large entanglement, as has been proven in refsP^EIl 
These models give again rise to long-range Hamiltonians, 
but they still very clearly demonstrate a strong distinction 
between correlations and entanglement. 

C. General gapped local spin models 

We will now turn to the discussion of general 1-D gapped 
spin models with local interactions, where each site supports 
a (i-dimensional quantum system. As it is stated rigorously 
in the theorem below for such models an area law always 
holds " ' . The proof is deeply rooted in the existence of Lieb- 
Robinson bounds which have also been essential in the proof 
of the expo nential decay of correlation functions in gapped 
local model^nMH.v 

As we allow for arbitrary d, it is sufficient to consider 
Hamiltonians on the chain L = {1, . . . , N} that have inter- 
actions only to nearest neighbors. Then 

where j+i is supported on sites j and We also impose 
a constraint of finite-interaction strength in that the operator 
norm i+i|| < J for some J > 0. Then re f.^^ finds: 

Theorem 5 (Area law for gapped spin chains) Consider a 
local Hamiltonian H as in eq. ( |2i| ) with finite interaction 
strength. Suppose H has a unique ground state with a spectral 
gap /S.E to the first excited state. Let us as before consider the 
block I = {!,..., rt}. Then, 

S{pi) < 5,„«, = coein(60 ln(d)26« (22) 

for some numerical constant cq > of order unity, and where 
^ = max(2v / AE , ^c)' v is the velocity of sound and (,c > 
is of the order of unity. 

The proof of this statement is quite intricat^~ -'and well be- 
yond the scope in this article. At its heart is the way locality 
enters by virtue of the Lieb-Robinson Theorem. It is a state- 
ment on the existence of a speed of sound in local Hamiltonian 
systems with finite-dimensional constituents: Let us imagine 



^ This result is compatible with an earlier result of an area law in 1-D gapped 
quantum field theories, based on the c-theorem presented in ref!^. This 
work also connected the role of the boundary points between regions / 
and O with the cluster decomposition in quantum field theory. In gapless 
system with open boundaries, the entropy is then half of the one in the 
situation of having periodic boundary conditions. 



we single out two disjoint sets X, Y from a lattice, and con- 
sider observables A and B that have support only on X and 
Y , respectively. Then [A, B] = 0. If we evolve A with time 
under a local Hamiltonian H it is no longer exactly true that 
A{t) and B commute: A{t) will be significantly supported on 
more and more sites, "melting away", and developing a long 
tail in support. For short times or large distances between 
sets X and Y, the commutator of A{t) and B will be very 
small. How small ex actly is governed by the Lieb-Robinson 
Theorer^^^^^^. 

Theorem 6 (Lieb-Robinson-Theorem) Let H be as in eq. 
( |27p a local Hamiltonian with a finite interaction strength. 
Then there exists a velocity of sound f > and fi,c> such 
that for any two operators A and B with support on disjoint 
sets X and Y we have that 

\\[A{t),B]\\<c\\A\\ \\B\\e^p{^^l{distiX,Y)-v\t\)), 

(23) 

where the distance between sets is taken to be dist{X, Y) = 
miiiigxjey (K ~ '^"'^ where 

A{t) = e'^*yle~'^*. 

The velocity v is of order J. 

This statement, natural as it may seem when viewed with 
a reasonable physical intuition, is a rigorous, and profound 
statement on how locality manifests itself in quantum lattice 
systems. From this bound, the deca y of co rrelation func- 
tions in gapped models can be proverPI^E^, an area law as 
above, U^, the quantization of the Hall conducance for inter- 
acting electrongiSI, as well as statements concerning propa- 
gation of quantum information and correlations through local 
dynamics-''.* Lieb-Robinson bounds also feature in the proof 
of a higher-dimensional Lieb-Schultz-Mattis theore m-'- 

We will later, in Subsection |IV.G| encounter another conse- 
quence of the Lieb-Robinson theorem, namely that quenched 
non-equilibrium systems generically satisfy area laws when 
starting from a product state and undergoing time evolution 
under a local Hamiltonian. This perspective receives a lot of 
attention in the context of non-equilibrium dynamics of quan- 
tum many-body systems. Here, the Lieb-Robinson is also the 
basis for the functioning of numerical light cone methods to 
study time evolution of quantum many-body systems Light- 
Cone,MPS,NonMPS,NonMPS2, in which effectively, only 
the essential part inside the causal cone is simulated. 

D. Results from conformal field theory 

In critical models the correlation length diverges and the 
models become scale invariant and allow for a description in 
terms of conformal field theories. According to the universal- 
ity hypothesis, the microscopic details become irrelevant for a 



The assumption that we have a spin system, meaning finite-dimensional 
local constituents, is crucial here. 
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number of key properties. These universal quantities then de- 
pend only on basic properties like the symmetry of the system, 
or the spatial dimension. Models from the same universality 
class are characterized by the same fixed point Hamiltonian 
under renormalization transformations, which is invariant un- 
der general rotations. Conformal field theory then describes 
such continuum models, which have the symmetry of the con- 
formal group (including translations, rotations, and scalings). 
The universality class is characterized by the central charge c, 
a quantity that roughly quantifies the "degrees of freedom of 
the theory". For free bosons c = 1, whereas the Ising univer- 
sality class has c = 1/2. 

Once a model is known to be described by a conformal 
field theory, powerful methods are available to compute uni- 
versal properties, and entanglement entropies (or even the full 
reduced spectra) of subsystems.^ This approach applies for 
1 + 1-dimensional systems, that is with one spatial dimension. 
In the seminal work ref."^^'the entanglement entropy in 1 + 1- 
dimensions has been calculated, see also refs.**^ The work 
refs.E^ put this into a more general context, and also allow 
for non-contiguous regio ns /. The local spectra of the reduc- 
tions are discussed in r efs j'^^ l '^^ J^. Block-block entanglement 
is also discussed inref.™ For a short non-technical review, 
see ref."^**. 

Starting point of the computations is the observation that 
powers of the reduced density matrix p" for any positive in- 
teger n can be computed. The series tr[p"] — Xj{pj)'^ 
is absolute convergent and analytic for all Re(n) > 1. The 
derivative exists, and hence one can make use of 

lini-^|-tr[p?] 
ri\,i In 2 an 

to compute the entanglement entropy. This procedure is typ- 
ically referred to as "replica trick". This leads in 1 + 1- 
dimensions to the expressiori^ 



S{pi)^-\og^{lla)+0{l), 



(24) 



where c is as above the central charge, I is the length of a 
single interval forming region /, a is an ultraviolet cutoff, cor- 
responding to a lattice spacing, to avoid an ultraviolet diver- 
gence, cp. eq. ( 10 1. The above constant C is hence nothing 



but c/3. This divergence i s also removed by using the mutual 
information"*^, see Section V.B The offset constant in eq. (24 1 



is non-universal. So the logarithmic divergence of the entan- 
glement entropy in the length of the interval is recovered here. 
From the expression given in ref.— for p], one also finds for 
the Renyi entropies for a > 1 



SaiPl) 



-(1 

6 ^ 



l/a)log2(Va)+0(l). 



' Conformal field theory provides — in this context specifically in 1 + 1- 
dimensions — a powerful repertory of methods to compute quantities that 
are otherwise inaccessible especially for non-integrable models. From a 
mathematical physics perspective, it is the lack of a rigorous proof of the 
relationship between the lattice model and the conformal field theory that 
makes such a treatment, pedantically speaking, non-rigorous. 



If one is close to the critical point, where the correlation length 
^ > is large but finite, one can often still effectively describe 
the system by a conformal field theory. One then obtains for 
the entanglement entropy^^ (compare also ref.-^^ 



E. Disordered spin chains 

Natural systems will generally exhibit a certain amount of 
quenched disorder which means that the model parameters are 
drawn randomly and the resulting correlation functions or en- 
tanglement entropies ]E[5(/9/)] have to be considered as being 
averaged over the a priori distributions, with average E. The 
critical behavior of quantum spin chains with "quenched" dis- 
order is remarkably different from its counterpart in the corre- 
sponding pure case, in several respects. Hence, it is only nat- 
ural to ask whether the scaling of the entanglement entropy is 
influenced by having some disorder in the model. This ques- 
tion has first been posed in ref.E^for the spin-1/2 random 
anti- ferromagnetic Heisenberg model. 



H 



E 



with { Jj} drawn from a suitable continuous distribution. The 
low energy properties of this model, along with the ran- 
dom XX model, are described by what is called a random- 
singlet phas^^. Using a real-space renormalization group 
approach'^^, the intuition can be developed that in this phase, 
singlets form in a random fashion, distributed over all length 
scales. The entanglement entropy of a sub-block is hence 
obtained by effectively counting the singlets that cross the 
boundary of the sub-block. This intuition is further developed 
in ref.G^. Within the framework of a real-space renormaliza- 
tion group approach — it is shown that the averaged entangle- 
ment entropy for a large class of disordered models scales hke 



0(1). 



In this class one hence observes universal behavior in the scal- 
ing of the averaged entanglement entropy. The intuition elabo- 
rated on above is further corroborated by work on the random 
anti-ferromagnetic XXZ chaiiV^, described by a Hamiltonian 



H 



3 "3+1 



AjCr|cr|+i) , 



where {Jj} are positive, uncorrected random variables 
drawn from some probability distribution and the uncorrected 
anisotropy parameters {A^} are also taken from a probability 
distribution. In this work, the observation is further explored 
that the scaling of the averaged entanglement entropy can be 
universal, even if correlation amplitudes are not, in that they 
would manifest themselves only in non-leading order terms in 
the entanglement entropy. This intuition is also further cor- 
roborated in refs! '^ ' i on the entanglement entropy in a 2-D 
situation. 
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From a fully rigorous perspective, the entanglement en- 
tropy in the random Ising model — for which ref "^'^ finds a 
scaling with the effective central charge of ln(-\/2) — has re- 
cently been revisited with methods and ideas of percolation 
theory^^^^. This approach is more limited than the Fisher- 
Hartwig techniques in terms of the class of models that can be 
considered — the Ising model only — but is more powerful in 
that also disordered systems with no translational invariance 
can be considered. 

Theorem 7 (Non-translationally invariant Ising model) 

Consider the Ising model 



H 



E 

j,keL,dist{j,k) = l 



where Xj ^ > and Sj > are the spin-coupling and external 
field intensities, respectively, which may depend on the lattice 
site in a non-translationally invariant system. The total num- 
ber of sites is N = 2m + n + 1, with {1, . . . , n} being the 
distinguished region. Then there exist 7, a, C with properties 
as in the subsequent footnote If^ > 4 In 2, then there exist 
constants ci , C2 > depending only on 7 such that 



S{pi) < ci log2(n) +C2, 



(26) 



for m > 0, so the entanglement entropy is at most logarithmi- 
cally divergent. 

The general picture that emerges is that the entanglement 
entropy scales as in the non-random case, but with a differ- 
ent prefactor in the logarithmic divergence. This seems natu- 
ral, as the disorder tends to "localize" excitations, and hence, 
with faster decaying correlations one would expect less en- 
tanglement to be present in the system. Yet, there are excep- 
tions: Cases in which one does find a logarithmic divergence, 
but with a larger prefactor compared to the non-random case. 
This includes the random quantum Potts model with spin di- 
mension d: Here, for the very large dimension of d > 41 one 
finds a larger factoi '^'. The exploration and complete classi- 
fication of the role of disorder to entanglement properties of 
ground states — including non-critical and higher-dimensional 
models — remains an interesting challenge. 



F. Matrix-product states 

Matrix-product states (MPS) play a very central role in the 
context of area laws for the entanglement entropy. They form 



'"Let A, (5 S (0,00) and write 9 = X/5. There exist constants o,C G 
(0, oo) — depending on 6 only — and a constant 7 = 7(6) satisfying < 



7 < 00 if 6 < 1, sucli tliat, for all n > 1, 

IIp^-PmII <min{2,Cn"e-^'^}, 



2 < Af < M. 



(25) 



Here denotes the reduced state of n sites in a system of total size A'^. 
One may find such a 7 satisfying 7 — ^ 00 as S \j 0. 



the class of states that is at the root of the workhorse of sim- 
ulating strongly correlated quantum chains — DMRG. This 
link will be elaborated upon in detail in Section VI Here, 



we focus on the entanglement and correlation properties of 
MPS. In the original sense, MPS are states defined on quan- 
tum chains consisting of N sites, each constituent being a d- 
level system. There are several ways of defining and intro- 
ducing MPS, the relationship of which may not be entirely 
obvious. This is also the reason that it was left unnoticed for 
some time that MPS — as being generated in DMRG — and 
finitely correlated stated — as being considered in the math- 
ematical physics literature — are up to translational invariance 
essentially the same objects. 

One way of looking at MPS is via a valence-bond picture: 
For each of the constituents one introduces a virtual substruc- 
ture consisting of two particles. Per site with Hilbert space 
<C^, one associates a Hilbert space (D^ (g) (D^ for some D. 
This D is sometimes referred to as the dimension of the cor- 
relation space, or D called the auxiliary or virtual dimension. 



Two virtual particles 
per physical site 




|'0)= J2 tr[AWH...ylWM]|ii,...,i^) 



il,...,tw = l 



FIG. 2 The valence-bond picture underlying matrix-product states 
as explained in the main text. 

These I?-dimensional virtual systems are thought to be pre- 
pared to a maximally entangled state with each one particle of 
each of the neighbors, arranged on a ring (see Fig.|2]). In other 
words, one starts from a pure state defined by the state vec- 
tor IV'd)®", where we have defined the maximally entangled 
state vector as 



(27) 



fe=i 



Then, one applies a local linear map to each of the pairs of sys- 
tems associated with every physical constituent in the center 
of the chain. 



(28) 



j = l a,b=l 



where k = 1, . . . ,N. This procedure will prepare a certain 
class of states: Indeed the MPS. We may conceive for each 
site k L the collection of complex numbers A^^l [j] as the 

elements of d matrices v4''''' [1], . . . , A^'^-' [d]. For a quantum 
spin chain with d = 2, we hence simply have two matrices 
A^'^) [1] , a'*') [2] per site. This procedure of locally projecting 
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to the physical dimension d gives rise to state vectors of the 
form 



(29) 



This is the most frequently used form of representing ma- 
trix product states. For open boundary conditions in a chain 
L ~ {1, . . . , N}, A'^^ [ii] and A'^^^ [ipf] are row and column 
vectors, respectively. MPS are described by only a number of 
parameters polynomial in the system size, 0{dND^), in con- 
trast to the scaling of the dimension of the full Hilbert space 
(C*)^^, which is exponential in TV." 

The particular construction of MPS immediately shows that 
MPS satisfy an area law. In fact, it follows trivially from their 
definition (see also fig. [2]i that 

S{pj)<2\og,{D), 

so the entanglement entropy is always bounded from above by 
a constant in N. MPS have hence an in-built area law prop- 
erty. As we have shown earlier the ground states of a variety 
of Hamiltonians exhibit exactly such an area-scaling when the 
system is non-critical but a logarithmic divergence when the 
system is critical. This already suggests that MPS may be a 
good description for ground states of non-critical systems but 
that this description may become less efficient in critical sys- 
tems. Indeed, it will be discussed and highlighted later in this 
article that systems satisfying an area law can be economi- 
cally represented as MPS so that MPS with a small auxiliary 
dimension D can indeed typically approximate ground states 
of local Hamiltonians. 



G. Single-copy entanglement 

The entanglement entropy — occupying center stage in this 
article — quantifies entanglement in a very precise sense: For 



" If one allows D to (exponentially) grow with the system size, one 
can easily show that actually every state vector from (C'')®^ can be 
represented as a MPS of the form as in eq. {29\ . It is important to note 
that MPS can not only be described with linearly many parameters in 
the system size: One can also efficiently compute local properties from 
them, which is a property not merely following from the small number of 
parameters to define them. For expectation values of observables having a 
non-trivial support on sites k, . . . , k + I £ L, we find {Sk ■ ■ ■ Sk+i) = 



(1) 



-,(fe-i)„(fe) 



..E 



(fc+0 „(fc+! + l) 



^1 



.(iV)l 



where the 



tt[E'^ ... ^1 ■■■ ^Sk+, 

transfer operators are defined as 

4''= E 015|fe>(AW[fc]®(A«[i])*), 

the star denoting complex conjugation. The decay of correlation functions 
can also be studied: If all matrices are the same per site, A^-^^ [j] =: A[j] 
for all j 6 L, and similarly define Ei , then one finds 

{SkSk+i) - {Sk){Sk+i) = 0{\X2{Ei)\'-^), 

where A2 {Ei ) denotes the second to largest eigenvalue of the transfer op- 
erator of the identity Ei . 



pure States it is the distillable entanglemen^^^^, so the rate 
with which one can locally extract maximally entangled pairs 
from a supply of identically prepared system. Specifically, 
local refers here to a subsystem / of the system, but to a col- 
lective operation on many identically prepared states. In a 
quantum many-body system, needless to say, this means that 
one performs operations that are local to all constituents in / 
collectively in all specimens at hand. 

When having the entanglement content in mind, one can 
equally reasonably ask how much entanglement is contained 
in a single spin chain. The concept of single-copy entangle- 
ment grasps this notion of distilling entanglement from a sin- 
gle specimen of a quantum spin chain with certainty. 

If D is the largest integer such that one can determinis- 
tically transform a state into the maximally entangled state 
|V'-d)(V'£i| (see eq. (27 1) by local operations and classical com- 



munication (LOCC), i.e. 



P ^ \^D){iw\ 



(30) 



one assigns the value Ei — log2(£') to the state as its single- 
copy entanglement. For pure states, such transform ations on 
the level of specimens are perfectly well understoocP^MHIllo] 
and is linked to the well-established theory of majorization in 
linear algebra'^"^. For our present purposes, for a pure state 
p = we find that eq. (30i holds if and only if ||p/|| < 

l/D. Hence, 

E,{pi)^\og^{W\pi\\-^\). 

This, in turn, means that the single-copy entanglement can be 
derived from the a-Renyi entropy of the reduction in the limit 
of large a. A quite surprising insight is that in critical systems, 
we do not only find a local spectrum leading to the logarith- 
mic divergence of the entanglement entropy. But that there is 
more structure to the spectrum, governing all of its Renyi en- 
tropies. For example, for quasi-free models, we find that once 
the entanglement entropy diverges, so does the single-copy 
entanglement, with a prefactor that is as ymptot ically exactly 
half the value of the entanglement entropjE^ill 

Theorem 8 (Single-copy entanglement) Consider a family 
of quasi-free fennionic Hamiltonians as in Theorem^ Then, 
whenever the entanglement entropy scales as 

^(p/) = Ciog2N + o(i), 

for some constant ^ > 0, then the scaling of the single-copy 
entanglement is found to be 

i?l(p/)-fl0g2(")+0(l). 

This means that exactly half the entanglement can be dis- 
tilled from a single critical chain than what is available as a 
rate in the asymptotic setting'^** '^''). This finding has also 
been corroborated by the behavior of all critical models for 
which the local spectra can be described b y th eir conformal 
field theory in quite some generalitySIl Ref.12^ studies Renyi 
entropies in boundary critical phenomena, and hence also ar- 
rives at a relationship between the entanglement entropy and 
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the single-copy entanglement. Ref. considers the entropy 
loss along the renormalization group trajectory driven by the 
mass term in free massive theories, and discusses also the 
single-copy entanglement in such situations. Ref.'^^ studies 
the situation of single-copy entanglement in the situation of 
bipartite systems between blocks when there is a gap of a fi- 
nite number of sites between the two blocks. Interestingly, 
there are critical models in which the single-copy entangle- 
ment still diverges in this sense. 

H. Summary of one-dimensional systems 

In a nutshell, the situation in one-dimensional translation 
invariant models is quite clear: If a system is local and gapped, 
an area law always holds rigorously. In many specific mod- 
els, prefactors can be computed. In contrast, if the interac- 
tions may be long-ranged area laws may be violated. For crit- 
ical lattice models for which one can directly evaluate the en- 
tanglement entropy, a logarithmic divergence is encountered. 
This picture is supported by the findings of conformal field 
theory. The situation will be less transparent and more intri- 
cate in higher-dimensional models. In any case, in the light 
of the previous findings, one may be tempted to formulate the 
following conjecture on the numerical bound on the right hand 
side of the previously discussed area law: 

Conjecture 1 (Area bound in one-dimension) There exists 
a function f : — >■ R+ — equipped with further suitable 
properties — such that in any gapped one- dimensional model, 
we have 

S{pi) < f{v/AE), 

where AE is the spectral gap and v is the speed of sound as 
used in the Lieb-Robinson bound. 

Indeed, while most explicit studies do indicate a behavior 
linear in ln(l/Ai?) of the entanglement entropy — the above 
mentioned quasi-free models — one can construct models'-*^ 
for which one finds a dependence which is polynomial in 

1/AE. 

IV. AREA LAWS IN HIGHER DIMENSIONS 



boundary of the region / is a non-trivial object. Now we are 
in a position to approach the question: Given a ground state of 
a quantum many-body system, does the entanglement entropy 
of a subregion / fulfill an area law? This question has been 
initiated in refs.'^^'^, where also a numerical answer has been 
found. 

The answer to this question for ground states is very much 
developed in case of — once again — quasi-free bosonic or 
fermionic models. Even in such systems, the rigorous an- 
swer to this question will turn out to be technically quite 
involved. The reason for these technicalities is essentially 
rooted in the very fact that one distinguishes a subregion, thus, 
e.g., breaking translational symmetry of translationally invari- 
ant systems, and analytical methods are hard to come by. The 
first rigorous higher-dimensional area law has been proven — 
to the knowledge of the authors — in refs.^ ^, w ith refinements 
for arbitrary harmonic interactions in ref.E^I^S that for such 
bosonic free models the problem can be considered solved in 
all generality forming a "laboratory" of what one should ex- 
pect in general systems. 

For critical fermionic model4^ ^^ I ^^ I ^^ I ^^^ J^, quite intrigu- 
ingly, one can find small violations of area laws: The area law 
is then only satisfied up to a logarithmic correction. In this 
section, we will discuss quasi-free models in great detail. Be- 
yond such quasi-free models, no rigorous results are known 
for states at zero temperature, with the exception of classes of 
states that satisfy an area law by their very construction, and a 
subsection will specifically be devoted to those. 

The models discussed here, however, do provide a clear 
intuition: Whenever one has a gapped and local model, and 
hence a length scale provided by the correlation length, one 
should reasonably expect an area law to hold. In cases where 
the number of eigenstates with vanishing energy-density is not 
exponential in the volume - a technical condition in its own 
right - one can even prove an area law with at most a logarith- 
mic correction from a sufficient decay of correlations'^^. The 
converse is not true, as we will see later, and one can have 
area laws even for critical systems in which the correlation 
length does not provide a length scale. For systems at non- 
zero temperatures, by contrast, the entropy of entanglement 
neither forms a meaningful measure of entanglement nor for 
quantum correlations. For appropriately defined measures of 
correlations, however, one can restore an area law which holds 
in generality for a large class of systems. 



For a chain, to satisfy an "area law" for the entanglement 
entropy means simply that it saturates with increasing block 
sizes. Needless to say, the notion of having entropic quan- 
tities scaling like the boundary area of a subregion becomes 
specifically relevant in case of higher dimensions: Then the 



Note that the single copy entanglement still grasps bipartite entanglement 
in a quantum many-body system. Identifying scaling laws for genuine 
multi-particl e entan glement is an interesting enterprise in its own right. 
Notably, refsEMIS) consider the geometric entanglement (the logarithm of 
the largest Hilbert-Schmidt scalar product with a pure product state) and 
relate it to the conformal charge of the underlying model. 



A. Quasi-free bosonic and fermionic models: Sufficient 
conditions for an area law 

We will follow the general description of ref.l^, where we 
think of the model being defined on a general lattice L speci- 
fied by a general simple graph. We consider quadratic bosonic 
Hamiltonians as in eq. Q and quadratic fermionic Hamiltoni- 
ans as in eq. (11 1. The key step is to relate correlation func- 
tions to entropic quantities. As before in the case of a har- 
monic chain, it is very involved to think of the entropy of en- 
tanglement itself. What comes to our help, however, is again 
that we can use the logarithmic negativity as an upper bound to 
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the entanglement entropy (see eq. (j6|l). The logarithmic nega- 
tivity is easier to treat analytically, as we can at all times refer 
to the full system, and not to subsystems /. In fact, we find 
that the logarithmic negativity can be bounded from above by 
the Li-norm of a submatrix of the covariance matrix^SI Por 
fermions in turn, the entropy may be bounded directly using 



the bound in eq. ( 20 1 



Theorem 9 (Entropic bounds from matrix norms) The en- 
tanglement entropy of ground states of quadratic bosonic 
Hamiltonians as in eq. Q satisfies 

s{pi)<En{pJ)<s\\t,\\ J2 Kp^pj)!- 

The entanglement entropy of unique ground states of 
quadratic fermionic Hamiltonians as in eq. Hjl satisfies 



s{pi)<2 J2 \ifU^) + {m\- 

iei.jeo 

This is a key tool towards proving the main theorem: We 
can reduce the evaluation of an entropic quantity to a counting 
argument over terms that can be evaluated from two-point cor- 
relators. Note that the use of the logarithmic negativity results 
in an important simplification of the problem. This shows that 
ideas from quantum information theory indeed help in find- 
ing proofs of statements of the scaling of the entanglement 
entropy. 

We are now in the position to state the bound of the scaling 
of the ground state entanglement in the bou ndary area s{I), 
eq. ([l]), of the distinguished region jEISIlII it is remarkable 
that merely the decay of two-point correlations matter here, 
and that even some critical models will give rise to an area law, 
as long as the algebraic decay of correlations is sufficiently 
strong. 

Theorem 10 (Quadratic Hamiltonians on general lattices) 

Let rj ^ V + 1 + 2e, e > 0, and assume that the ground state 
is unique and fulfills for i, j £ L, i ^ j, and some constant 
Ko 



dist'^{i, j) 



> 



\{PtP3)\ 



for bosons, 
for fermions. 



Then 



S (pi) <KoCvCil +£)■<!) X 



Ir^;!! for bosons, 
1 for fermions, 



where ^ is the Riemann zeta function and the constant c-p de- 
pends only on the dimension of the lattice. 

A general version of what one should expect to be true pro- 
vides the connection to the spectral gap: For gapped models 
the correlation functions decay exponentially with the graph 
theoretical distance. One cannot apply the Lieb-Robinson 
Theorem to prove this, unfortunately, as the involved opera- 
tors are unbounded. Hence, a technique that is applicable to 
describe clustering of correlations in such models had to be 



developed. The ideas of the proof go back to ref.El, general- 
ized to arbitrary lattices in ref .Emu. Key ideas of the proof are 
polynomial approximations in the sense of Bernstein's Theo- 
rem. For a thorough discussion of clustering of correlations 
in translation-invariant harmonic systems, see ref.^*". For 
general lattices and gapped quadratic bosonic and fermionic 
Hamiltonians, one finds that two-point correlation functions 
decay exponentially. Together with the above theorem this 
leads to an area law whenever the model is gapped: 

Corollary 1 (Area law for gapped quasi-free models) The 

entanglement entropy of ground states of local gapped models 
of the type of eq. ^ for bosons and of eq. ( |77| ) for fermions 
for arbitrary lattices G — (L, E) and arbitrary regions I 
satisfies for a suitable constant ^ > 

S{pi) < is{I). 



B. Logarithmic correction to an area law: Critical 
fermions 

What can we say about situations in which the previous suf- 
ficient conditions are not satisfied? Specifically, how is the 
scaling of the entanglement entropy modified in case of criti- 
cal fermionic models? This is the question that will featur e in 
this subsection. Following the bosonic result in refs.l^^EllI, the 
entanglement entropy in fermionic models was first studied 
in ref.^^^^ for cubic lattices. Here, the quadratic bound in eq. 
( 20 1 plays an important role, to relate bounds to the entropy 
to feasible expressions of the covariance matrix of the ground 
state. Here, not quite an area law, but only one up to a loga- 
rithmic correction is found. The results can be summarized as 
follows: 



Theorem 11 (Violation of area laws for critical fermions) 

For a cubic sublattice I = {l,...,n}^^ and an isotropic 
quasi-free model as in eq. with a Fermi sea of non-zero 
measure and a finite non-zero surface there exist constants 
Co, Ci > such that the ground state fulfils 



Con 



ln(n) < S{pi) < Cin^-Mn^(?i) 



The stated lower bound makes use of the assumption that the 
Fermi surface is finite (and of a technical assumption that the 
sets representing the states cannot have nontrivial irrelevant 
directions); assumptions both of which can be removed^''. 

This fermionic quasi-f ree case akeady exhibits a quite com- 
plex phase diagranP^EHl. At the same time, ref.^* formulated a 
similar result, based on a conjecture on the validity of Fisher- 
Hartig-type scaling result for higher dimensional equivalents 
of Toeplitz matrices, as was further numerically corroborated 
in ref.'^. A logarithmic divergence is not directly inconsistent 
with the picture suggested in a conformal field theory setting, 
as relativistic conformal field theories do not have a Fermi 
surfacd'^-5J. It is still intriguing that critical fermions do not 
satisfy an area law, but have logarithmic corrections. In this 
sense, critical fermionic models could be said to be "more 
strongly entangled" than critical bosonic models. 
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C. Difference between critical fermions and bosons: IHalf 
spaces 

The scaling of block entropies for bosons and fermions 
in higher spatial dimensions hence exhibit remarkable differ- 
ences. Let us consider the case of a cubic lattice of nP sites 
with periodic boundary conditions and / = {1, . . . ,m} x 
{1, . . . , n}^^^^ (w.l.o.g. we distinguish the first spatial di- 
mension). Then one may transform the Hamiltonian to a sys- 
tem of mutually uncoupled one-dimensional chains using a 
unitary discrete Fourier transform. After this decoupling pro- 
cedure the entanglement between / and O is given by a sum 
of the entanglement between the sites /' = {1, . . . , to} and 
O' — {to + 1, . . . , n} of the ttP^^ individual chains 



Sjpi) 
s{I) 



,-D-l 



1 



n 

■E 

1=1 



S{p},). 



We start with a discussion of fermions and focus on the 
isotropic setting (B = in eq. ([TT])). After taking the limit 
n — ?► oo, the asymptotic behavior in to of the entanglement 
S{pj, ) can then be read off Theorem to yield the following 
statement (for technical details see refP^: 

Theorem 12 (Prefactor for fermionic half spaces) 

Asymptotically, the entanglement entropy of fermionic 
isotropic models of half spaces satisfies 



lim 



S{pi) ^ l0g2("x) 

s{I) ( 



s=l 



0(1). 



Here, 



_ |{0e [0,2^)^-1 1 

(27r)^-i 

is the integral over individual chains (j) with s discontinuities 
a{(j)) in their symbol. 

Hence, one encounters a logarithmic divergence in to of the 
entanglement entropy and the pre-factor depends on the topol- 
ogy of the Fermi-surface: The symbols exhibit discontinuities 
on the Fermi-surface. If the Fermi surface is of measure zero 
(i.e., the set of solutions to = 0, (/) G [0, 27r)^, is count- 
able, as, e.g., in the critical bosonic case discussed below), we 
have Vs = and the system obeys the area law. 

For the bosons, we discuss the important case of to = n/2 
for the Klein-Gordon Hamiltonian as in eq. (j9]). After the 
transformation to uncoupled chains, one finds Hamiltonians 
for the individual chains that correspond to a nearest-neighbor 
coupling matrix X of the form as in Theorem [T[ which yields 



En{p,I) 
,s(/) 



dip 



in the limit n ^ oo. This expression is independent of the 
mass TO and finite: For 2? = 2, it evaluates to logj (3+2\/2) /4 
and similarly for V > 2. Hence, despite being critical, the 



system obeys an area law, in contrast to the fermionic case 
(for TO — n/2 the entanglement for a critical fermionic system 
would diverge in n). 

Hence, in quasi-free critical models, it matters whether a 
system is bosonic or fermionic when it comes to the question 
whether or not an area law holds. The above results confirm 
the numerical analysis of ref.'^"^' for critical bosonic theories, 
and of ref.G^ for two-dimensional fermionic systems. Mo- 
tivated by these findings, numerically studies the non- 
leading order terms of an area law in nodal fermionic systems: 
It is found that in non-critical regimes, the leading subarea 
term is a negative constant, whereas in critical models one 
encounters a logarithmic additive term. A lesson from these 
higher-dimensional considerations is that the simple relation- 
ship between criticality and a violation of an area law is hence 
no longer valid for local lattice models mV > 1. 



D. Entanglement in bosonic thermal states 

In this subsection, we briefly discuss area laws for notions 
of entanglement in Gibbs states. 



tr[exp(-/3i/)] 



for some inverse temperature /? > 0. The second moments 
matrix, the covariance matrix, is then found to be F = © 

r, = x-i/2(xi/2pxi/2)'/'(i + G)x-'/\ 

G - 2(cxp[/3(xi/2pxi/2)'/'] - iy\ 

Using the the methods of refsl ^' l ^^ ^one again finds the suitable 
decay of correlations, which can be translated into an area law 
for the e ntanglement content. Here, the result — taken from 
-is stated in terms of the logarithmic negativity. 

Theorem 13 (Entanglement in thermal bosonic states) 

The logarithmic negativity of thermal states of quadratic 
finite-ranged bosonic Hamiltonians as in eq. Q for 
[X, P] — satisfies En{p, I) < ^s{I)for a suitable constant 
^>0. 

Since the logarithmic negativity is an upper bound to 
the entanglement of formation and hence the distillable 
entanglement^^^ this implies an area law for these quanti- 
ties as well. It is important to stress that the entropy of a sub- 
region as such no longer reasonably quantifies entanglement 
between that subregion and the rest of the lattice: Even classi- 
cally correlated separable states will in general have a positive 
entropy of the reduced state. The latter quantity is then indeed 
extensive and fulfills a volume law, unlike the entanglement 
content. Area laws in thermal states have further been studied 
in detail in ref.^, where an emphasis has been put on identi- 
fying regions where the states become separable. Refs.l^SlSIl 
investigate thermal bound entanglement — entanglement that 
is not distillable — in bosonic quadratic and spin systems. 
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E. Results from conformal field theory 



resulting state vector as such hence becomes 



In systems with more than one spatial dimension, the sit- 
uation is more intricate, and there is no general expression 
known for entanglement entropies in d + 1-dimensional con- 
formal field theories. For interesting steps into a description of 
systems with two spatial dimensions in the framework of con- 
formal field theory see refs.^" For a class of critical mod- 
els in two spatial dimensions (including the quantum dimer 
model), it is found that S{pi) — 2fs{L/a) + cg\n{L/a) + 
0(1), where L is the length of the boundary area, fg is an 
area law prefactor that is interpreted as a boundary free en- 
ergy, and g is a coefficient that depends on the geometric prop- 
erties of the partition into / and O. That is, in addition to a 
non-universal area law, one finds a universal logarithmically 
divergent correction. For a further discussion of steps towards 
a fuU theory of entanglement en tropies in c? + 1-dimensional 
conformal field theories, see refs.l^SESl 



F. States satisfying area laws by construction: Projected 
entangled pair states, graph states, and entanglement 
renormalization 

In this section, we will discuss classes of states that have 
the area law already built into their very construction. In 
this sense, they grasp the entanglement structure of local 
higher-dimensional models. These are projected entangled 
pair states, so matrix-product states in higher dimensions, and 
states from entanglement renormalization. They are designed 
to be variational states well-approximating true ground states 
of local many-body systems: As was already true for matrix 
product states, they form a complete set of variational states. 
Yet, typically, for a much smaller, polynomial or constant, 
number of variational parameters, they often deliver a very 
good approximation. In projected entangled pair states, local- 
ity is respected in just the same way as for MPS. Entangle- 
ment renormalization, in turn, is based on a scale-invariant 
tree structure, intercepted by disentangling steps, which in 
higher dimensions nevertheless leads to an area law for the 
entanglement entropy. 

Projected entangled pair states (PEPS) can be thought of 
as being prepared as MPS in higher dimensional cubic lat- 
tices L = {1, . . . , N}^^, or in fact to any lattice defined 
by any undirected simple graph G ~ {L,E). In this va- 
lence bond construction, one again associates a physical space 
with Hilbert space (D'' with each of the vertices L of G. 
Then, one places a maximally entangled pair of dimension 
D X D (see eq. ( 27 1) for some positive integer D between 
any two vertices that are connected by an edge e £ E. For 



E 

ll,l2,...,4lL|=l 



,«|L|) 



a cubic lattice, one hence starts from a cubic grid of maxi- 
mally entangled state vectors. Then, one applies a linear map 
p(fc) : <C^ (g) ■ ■ ■ (g) (C^ ^ (C^ to each physical site, as 

j = l ll,...,i|Si(fc) 

Here, 1 51(^)1 is the vertex degree of the vertex k E L. The 



where C denotes a contraction of all higher-order tensors with 
respect to the edges E of the graph. This amounts to a summa- 
tion over all indices associated with connected vertices. The 
objects A'^'^^ are hence tensors of an order that corresponds 
to the vertex degree of the lattice (a second order tensor — a 
matrix — for a one-dimensional chain, a three order tensor in 
hexagonal lattices, a fourth order tensor in cubic lattices with 
2? = 2, and so on). This construction is the natural equivalent 
of the valence bond construction for matrix-product states as 
explained in eq. (28 1. This ansatz as such is the one of ten- 



sor product states that is due to ref.'^^ which in turn is gen- 
eralizing earlier work on AKL T-type valence bond states in 
two dimensions in refs.E^3S3 xhe generated class of states 
is referred to as projected entangled pair state s^^^ , reflect- 
ing the preparation procedure. PEPS states are sometimes 
also in higher dimensions simply referred to as matrix prod- 
uct states^^^. This ansatz has proven to provide a powerful 
and rich class of states. Importantly, ref.^'^ provides a first 
simulation method based on PEPS.'^ 

This class of states is complete, in that any state of a given 
finite lattice can be arbitrarily well approximated by such a 
state if D is sufficiently large. Clearly, to compute local ob- 
servables in such an ansatz, one has to contract this instance 
of a tensor network which in 2-D is actually computationally 
hard.'"* It is however possible to provide approximation tech- 
niques, related to the DMRG approach, that allow for the con- 
traction of the tensor network and then for the computa tion of 
the expectation values of local observablei^^iE^SlIIlEl] 

A particularly simple yet important subset of the projected 
enta ngled pair states is constituted by the so-called gr aph 
sfafejnmnillll They are instances of stabilizer states^^ 
which can be thought of as being prepared in the follow- 
ing fashion: On any graph G — {L,E), one associates 
each vertex with a (D^-spin. This spin is prepared in |+) = 
(|0) + |l))/2^/^. Then, one applies a p/ioie gafe 

U = |0,0)(0,0| + |0, 1)(0,1| + |1,0)(1,0| - |1, 

to each pair of vertices that are connected by an edge. This 
phase gate corresponds to an Ising interaction. Clearly, this 
construction makes sense for any simple graph, and this is 
a subset of the above projected entangled pair states. Graph 



" Note that 1-D MPS ba.sed on a suitable order of the constituents do not 
form a good approximation for 2-D models. This is essentially rooted in 
the observation that one should expect an area law for the entanglement 
entropy in gapped 2-D models. For an important early discussion of spectra 
of subsystems in 2-D integrable models see ref.^" and, for a more recent 
discussion of the implication on DMRG, ref.''^. 

In fact, it is known that the exact contraction of s uch a tensor network is 
contained in the complexity class ^P-completi^Hll Clearly, this means that 
no algorithm is known with polynomial running time. 
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states readily satisfy an area law by constructiorPS^E^HI as one 
merely needs to count the edges over the boundary of a dis- 
tinguished region to obtains the entanglement entropy, then 
obviously linear in the boundary area.'^ 

Gr aph sta tes may be generalized to weighted graph 
sfflfei^liiSUMI where the edges may carry a different weight, 
and in turn generalize to the ansatz of a renormalization al- 
gorithm with graph enhancement (RAGE), being a strict 
superset of matrix-product states and weighted graph states, 
one that can nevertheless be efficiently contracted. As the 
graph defining the (weighted)-graph state does not need to 
have the same structure as the graph of the physical system 
whose quantum state we would like to describe, (weighted)- 
graph states may describe volume scaling on the level of the 
physical system. This makes them particularly suitable for 
simulation of time evolution, where no area law can be ex- 
pected to hold. 

Yet a different class of many-body states with applica- 
tions in the simulation of quantum spin systems is given 
by the states generated by entanglement renormalization 
(MERA)^^^. This is a class of states the construction of which 
is inspired by a renormalization scheme. Consider a tree ten- 
sor network with the physical sites at the end. This can be 
efficiently contracted. Yet, when decimating, say, two spins 
of one layer to a single "superspin" in the next layer in a sin- 
gle step of a renormalization procedure , one loses information 
information about the state. The idea of a MERA ansatz is to 
allow for disentangling unitaries, effectively removing entan- 
glement from a state, before doing a renormalization step. 

More specifically, consider a cubic lattice L — 
{1, . . . ,iV}^^ in some dimension V, embodying sites. 
Each site j E L is associated with a physical system with 
Hilbert space (C*. The MERA is essentially a unitary ten- 
sor network of depth 0(hi(iV)), preparing \tp) from |0)®^. It 
consists of layers of isometrics — ^performing the renormaliza- 
tion step — and disentanglers, which minimize the entangle- 
ment in each step before the next renormalization step. This 
renormalization step may be labeled with a fictitous time pa- 
rameter. Each of the unitary disentanglers U G U{d"'^) in the 



The other side of the coin of the difficulty of actually contracting tensor 
networks, even if they correspond to states that approximate ground states 
satisfying area laws well, is that such states can have computational power 
for quantum computing. Indeed, certain graph states or cluster states — 
as they are called for a cubic lattice — are universal resources for quantum 
computing: Quantum computing can be done by merely applying local 
measurements onto single sites of such a cluster states, without the need of 
additional unitary control. This computational model — known as one-way 
computint'^^ — can also be understood as a teleportation scheme in virtual 
qubits^i^. The tensor networks that occur when performing Pauli measure- 
ments can still be efficiently contracted, but not under arbitrary measure- 
ments, leading to universal computation. The program of using general 
projected entangled pair states in quantum computing based on measure- 
ments only has been pursued in refsLMJSl giving rise to a wealth of new 
measurement-based quantum computational models. This also highlights 
how the disadvantage of having no classical efficient description can be 
made an advantage: One can at each instance of the computation not effi- 
ciently compute the outcome, but on a physical system realizing this model 
one could efficiently simulate any quantum computer. 



disentangling layer has a finite support on m sites. In the sim- 
plest possible realization of a MERA this would be m = 2. 
The unitaries can be taken to be different in each layer, and 
also different from each other within the layer Unlike PEPS, 
they do not give rise to strictly translationally invariant states, 
even if all unitaries are taken to be identical in each layer. 

Such a procedure can be defined for cubic lattices of any 
dimension T). In P = 1, one does in fact not observe an 
area law, but typically a logarithmic divergence of the entan- 
glement entropy, quite like in critical spin systems. Indeed, 
the MERA ansatz as a scale invariant ansatz is expected to 
be suitable to approximate critical systems well, and numer- 
ical sim ulations b ased on the MERA ansatz corroborate this 
intuitioiPMH. 

A precise connection between homoge- 
neous instances of a MERA ansatz and conformal field the- 
ory is established in refs.'^^ In more than one dimension, 
2? > 1, MERA satisfy again an area law, as a moment of 
thought reveals: On encounters asymptotically linearly many 
unitaries over a boundary that have entangling power, ren- 
dering the computation of an upper bound to the entangle- 
ment entropy a combinatorial problem. Despite this observa- 
tion, first numerical work on fermionic instances of a MERA 
ansatz appear to deliver promising results for strongly corre- 
lated systems''*^' ^2'**'. 

Theorem 14 (Area laws for PEPS, graph states, and MERA) 

For any finite dimension D of the virtual systems, the entan- 
glement entropy of a projected entangled pair state satisfies 
S{pi) < s{I)D, where as before s{I) denotes the surface 
area of I on a graph. Hence, also graph states with a fixed 
vertex degree satisfy area laws. A family of states from 
entanglement renormalization will also satisfies an area law 
for cubic lattices with V > 2, and a logarithmic divergence 
in P = 1. 

Interestingly, based on a PEPS description, one can con- 
struct critical models that still satisfy an area law in 2? = d^Sl 
resembling the situation for critical quasi-free bosonic sys- 
tems. The validity of an area law follows trivially from con- 
struction, so the technical part in the argument amounts to 
showing that a model is critical. In ref.'^'^'this is shown by em- 
ploying a quantum-classical-correspondence: Take a classi- 
cal two-body spin Hamiltonian of the form iJ(cri, . . . , ct^v) — 
Sdist(i j)=i h{ai,aj), ai — 1, . . . , d. This Hamiltonian will 
have at some inverse temperature /3 > a partition function 
Z = '^^e^^^^"''- From this classical partition function, 
a quantum state can be constructed by using the Boltzmann 
weights as superposition coefficients, 

(Tl,...,CTiV 

This state vector has the properties that for diagonal observ- 
ables, it gives rise to the same expectation values and corre- 
lation functions as the corresponding classical model does, it 
has a simple representation as a PEPS for D = d, and it is — 
as any PEPS — the ground state of a local Hamiltonian. The 
classical model can then be chosen such that the appropriate 
decay of correlation functions follows. This construction de- 
livers critical spin models that nevertheless satisfy an area law. 
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G. Quenches and non-equilibrium dynamics 

A physical setting that receives a lot of attention in 
the recent literature is the one of non-equilibrium dy- 
namics of quantum many-body systems. A specif- 
ically interesting setting is the one of a sudden 
qucncl J'^'^'^^''^^'^^'^^'^'*'^^™^^'^^'''"^'''*^'^^^'^^'*'^"^'^'^^^ Here 
the initial condition is the non-degenerate ground state of 
some local Hamiltonian H, with state vector Then, 
one suddenly (locally) alters the system parameters to a new 
Hamiltonian V. Since j?/') will typically no longer be an 
eigenvalue of H, one arrives at a non-equilibrium situation: 
The state vector's time evolution is simply given by 

Studies of instances of such complex non-equilibrium many- 
body dynamics and questions of the dynamics of quantum 
phase transitions are enjoying a renaissance recently, not the 
least due to the advent of the high degree of control over quan- 
tum lattice systems with cold atoms in optical lattices.'^ 

For finite times, infinite quench ed syste ms satisfy an area 
law in the entanglement entropjEEHSl (strictly speaking 
whenever one considers time evolution under local finite- 
dimensional Hamiltonians starting from product states). For 
finite systems this holds true for times that are sufficiently 
small compared to the system size over the speed of sound. 
The intuition is that when suddenly switching to a new Hamil- 
tonian, local excitations will be created. These excitations will 
propagate through the lattice, but — except from an exponen- 
tially suppressed tai l — at mo st with the Lieb-Robinson veloc- 
ity of Theorem |e[i vUHni] rphis IS yet agam a consequence 
for the approximate locality in quantum lattice systems, re- 
minding of the situation in relativity and implies that correla- 
tions can only slowly build up, resulting in an area theorem. 
In turn, such a quench does in general give rise to a linear 
increase in the entanglement entropy, a statement that is prov- 
ably correct, and has been encountered in numero us numerical 
studies on quenched non-equilibrium systems'^^^IIIIIllIIlIIoIl 
In fact, finite subsystems can locally relax in time, to appear 
as if they were in a thermal statJ^. These results may be sum- 
marized in the following statement. 

Theorem 15 (Area laws in non-equilibrium systems) Let 

\^) be a product initial state vector, and H a local Hamilto- 
nian. Then, for any time t > there exist constants Cq , ci > 
such that for any subset I the entanglement entropy of the 
time evolved reduction pj of pit) — e~'*^|'0)(^|e'*^ satisfies 

S{piit))<cos{I) + ci. (31) 

Specifically, this is true for any local Hamiltonian on a cu- 
bic lattice in dimension V. This means that for any constant 
time, the entanglement entropy satisfies what is called an area 



The interesting situation of locally perturbing the state and henc e gener- 
ating a non-equilibrium situation has also been considered in 
where an area law is always expected to hold. 



law. In turn, there are product initial state vectors \'ip) of 
one-dimensional spin chains, local Hamiltonians H, and con- 
stants C2, C3, C4, Lq, So, to > such that 

S{pi{t)) >C2t + C3, 

for L > Lq and s > Sq and < t < 04^3, for / = {1, . . . , s}. 

That is, for any fixed time t, one encounters an area law for 
the entanglement entropy, but the prefactor can grow linearly 
in time. In fact, by a suitable choice of blocks, one can show 
that a lower bound grows linearly in time! This fact is respon- 
sible for the hardness of simulating time evolution of quantum 
many-body systems using instances of the DMRG approach: 
to represent such states faithfully, exponential resources are 
then required. Similar bounds give rise to statements on the 
minimal time needed in order to prepare states with topologi- 
cal order using local Hamiltoniang^il 

There is an interesting localization effect of entanglement 
under quenched disorder, linking to the previous discussion 
on ground state entanglement in disordered systems. Whereas 
one obtains from Lieb-Robinson bound the estimate in time 

S{pi{t))<CQ\t\+ci 

for suitable constants co,ci > 0, in the disordered one- 
dimensional XY spin chain this bound is replaced by the 
tighter bound 

S{pi{t)) < co\n{N\t\) + ci, 

again for appropriate constantJ^H. This means that due to 
quenched disorder, the growth of entanglement is merely log- 
arithmic in time, not linear There is an intuitive explanation 
for this: The linear sound cone provided by the Lieb-Robinson 
bounds is replaced by a logarithmically growing or even a 
constant one, leading to a suppressed entanglement propaga- 
tion. A similar behavior is observed under time-dependent 
fluctuating disordei'2^. 



H. Topological entanglement entropy 

The topological entanglement entropy is a quantity that is 
constructed in a fashion that enables it to characterize quan- 
tum many-body states that ex hibit topological order, a con- 
cept introduced in refs F^^ I ^^ ^ (see also refs.'^^oilSS) On both 
sides of a critical point in a system undergoing a quantum 
phase transition, the quantum many-body system may have 
a different kind of quantum order; but this order is not nec- 
essarily one that is characterized by a local order parameter: 
In systems of, say, two spatial dimensions, topological order 
may occur. Topological order manifests itself in a degeneracy 
of the ground state manifold that depends on the topology of 
the entire system and the quasi-particle excitations then show 
an exotic type of anyonic quasi-particle statistics. These are 
features that make topologically ordered systems interesting 
for quantum computation, when exactly this degeneracy can 
be exploited in order to achieve a quantum memory robust 
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against local fluctuations. They even allow in theory for ro- 
bust instances of quantum computation, then referred to as 
topological quantum computatior^^^. 

The topological entanglement entropy is now designed 
as an instrument to de tect such topological order. In- 
troduce d in refs! ''^'' l ^"^^ i it received significant attention 
recentl>CISIIMilll32lIlim]_ The details of the relationship be- 
tween positive topological entanglement entropy and topolog- 
ical quantum order are discussed in ref.li^. 

In ref.''*** a disc in the plane / is considered with boundary 
length L. This disk is thought to be much larger than the cor- 
relation length, and it is hence assumed that an "area law" in 
the above sense holds. The entanglement entropy of pj will 
then have the form 



=aL-7 + 0(l), 



(32) 



where the last term vanishes in the limit L — > oo. The pref- 
actor a is non-universal and ultraviolet divergent. However, 
7 > is an additive constant which is universal and character- 
izes a global feature of the entanglement in the ground state. 
This quant ity i s referred to as topological entanglement en- 
tropy in ref.'"*". To avoid ambiguities when distinguishing the 
constant term from the linear one in eq. ( 32 1, ref. '^^ makes use 



of the following construction: The plane is divided into four 
regions, each of them being large compared to the correlation 
length. A, B and C are arranged as neighboring each other 
in three identical subparts of a disk. D is the exterior of the 
disk. The respective reductions to the parts are denoted as pA 
and pab to regions A and jointly A and B, respectively. The 
topological entropy S'xopo is then defined as 

5'topo = S{pa) + S{pb) + S{pc) (33) 
- S{pab) - S{pbc) - S{pAc) + S{pABc)- 

This is a linear combination of entropies of reductions, con- 
structed specifically in a way such that the dependencies on 
the length of the respective boundaries of regions cancel. It 
is not directly meant as an information theoretical quantity, 
although the differences of entropies resembling a mutual in- 
formation expression. Also, slightly different definitions with 
similar properties are conceivable, and indeed, the indepen- 
dent proposal of ref;^ makes use of an alternative combina- 
tion of entropies. The important aspect here is the above men- 
tioned cancellation of the boundary term. Taking the behavior 
as in eq. ( [32| for granted, one indeed finds 

•Sjopo = — 7- 

From the way S'xopo is constructed it is a topological invari- 
ant, and depends only on a universal quantity (unaltered under 
smooth deformations, as long as one stays away from critical 
points), and on the fashion how the regions are located with 
respect to each other, but not on their specific geometry (again 
assuming that the correlation length is much smaller the re- 
gions and does not matter). Interestingly, topological order is 
hence a global property that is detected by the entanglement 
entropy. This construction can also readily be used in numer- 
ical studies. The explicit computation how the entanglement 



entropy detects the presence of topological order in an actu- 
ally time-dependent model undergoing a quantum phase tran- 
sition from a spin-polarized to a topologically ordered phase 
has been systematically explored in ref.^^, further strength- 
ening the findings of ref.Q^. 

Since its proposal, this and related quantities have been 
considered in a number of contexts. A natural candidate to 
explore this concept is the toric code state of ref.''*': Con- 
sider for this a square lattice I — {1, ... ,n}^^ with periodic 
boundary conditions, and place the physical two-dimensional 
quantum spins on the vertices of this lattice.'^ This lattice is 
tiled into two sublattices of different color, red and white. Ev- 
ery white p and red plaquette s is then associated with one of 
the commuting operators 



A. 



jeds 



Bn 



n 



(34) 
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respectively, with non-trivial support on four spins each, 
where as before af , af denote the Pauli operators sup- 
ported on i. The Hamiltonian of the system — a local 
Hamiltonian — is then taken to be 



H = -Y.A,-Y.B, 



This is a gapped and frustration-free Hamiltonian. It is also 
straightforward to verify that for any closed path g the oper- 



ator rijgg '^j commutes with all operators in eq. (34i. The 
ground state manifold depends on the topology of me lattice 
and is in the chosen case four-fold degenerate. The topolog- 
ical entanglement entropy, evaluated for this toric code state, 
gives 7 = ln(2). The ground states can readily be cast into a 
PEPS language, as has been done in ref.l^. An analysis how 
topological order can be grasped in a language of entangle- 
ment renormalization or MERA has been performed in ref.El 
Indeed, the topological degrees of freedom can then be dis- 
tilled to the top of the tensor network. 

An equally important explicit and closely related model is 
the loop model on a honeycomb lattice of ref.l^. Ground 
states of more general string-net lattice models can also 
often be express ed in terms of remarkably simple tensor 
network^^DSlIMI Entanglement entropies of topological 
color code^ have been studied in ref. ' - J Equivalents of 
the topological entanglement entropy for finite temperature — 
where the very robustness c an be probed — have been consid- 
ered and introduced in refs.l^^nSl Notably, for Gibbs states 
it still makes sense to consider quantities of the type as in 



eq. (33 1, only with the respective entropies being replaced by 
mutual informations grasping correlations instead of entan- 
glement, as discussed in detail in Section [V] It is found that 
the interplay between thermal effects, topological order and 
the size of the lattice indeed give rise to well-defined scaling 
relations. 



Equivalently, one can place the physical spins on the edges and formulate 
the operators {As} and {-Bp} as being non-trivially supported on the re- 
spective four spins associated with vertices and plaquettes. 
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The study of entanglement entropies in fractional quan- 
tum Hall states in a spherical geometry has been initiated in 
j.gfti08j j-gfp43j Abelian Laughlin states as well as Moore- 
Read states have been considered, where also rigorous upper 
bounds for particle entanglement entropies have been derived. 
Particle partitioning entanglement in itinerant many-particle 
systems has been studied in ref.'^**'. The MPS representation 
of the Laughlin wave function has been derived in ref. The 
tolopogical entanglement of integer quantum Hall states has 
been computed in ref.^^. Topological entanglement Renyi 
entropies have been considered in ref.^^. Similar quantities in 
Chern-Simons theories — the best understood topological field 
theories — have been identified in ref.l^. The suggestion that 
the full spectrum of H in pj ~ e^^ should be considered 
to detect topological order has been proposed in ref.^SH. As 
being certified by this list of recent developments, studies of 
entanglement entropies as indicators of topological order are 
still under rapid development. 



from a common cause — the general locality of interactions) is 
still subject to debate.'^ 

This relationship has even been employed to take steps in 
computing the entanglement entropy in higher-dimensional 
conformal field theories: The AdS/CFT correspondence — 
relating a d + 2-dimensional anti de Sitter (AdS ) spa ce to 
a d + 1-dimensional conformal field theory (CFT)^^^^ — has 
been made us e of to study the Bekenstein formula in the AdS 
contexfSsnill see also ref.^^. In this way, the above formula 
is used as a tool to compute the geometric entropy in a plau- 
sible fashion in situations where the exact computation is not 
known to be possible using the tools of conformal field theory. 

The ho lographic principle — dating back to work in 
refs.GSHSll — goes even further, and suggests that generally, all 
information that is contained in a volume of space can be rep- 
resented by information that resides on the boundary of that 
region. For an extensive review, see ref.l^Sl. 



I. Relationship to blacl< liole entropy 

As mentioned before, one of the particularly intriguing mo- 
tivations for the study of area laws of the entanglement en- 
tropy is the suspected relationship to the area-dependence 
of the blac k hole entropy. The Bekenstein-Hawking area 
l^-^^Mnnni suggests that black hole carries an entropy that is 
proportional to its horizon area A, 

^ _ kc^A 

Hence, according to this relationship, the (thermodynamical) 
entropy of a black hole is just a quarter of its area measured in 
Planck units^'', i.e., when k^c — G — h^l. For the sum of 
this black hole entropy and the matter entropy 5Mattei a second 
law of thermodynamics is proposed to hold. Such a general- 
ized second law of thermodynamics led to the suggestion that 
one would have a "spherical entropy bound" for matter: In 
asymptotically flat spacetime, any weakly gravitating matter 
system would satisfy S'Matter < 2TrkEr/{hc), interestingly not 
containing G. E denotes the total mass energy of the system, 
whereas r stands for the smallest radius of a sphere that con- 
tains the matter system at hand. The range in which one can 
expect the validity of such a law is discussed in ref.^^. 

The linear relationship between the boundary area and 
the (thermodynamical) entropy — formally, the two equations 
look identical — suggests that one may expect a close rela- 
tionship between these area laws: On the on hand, for the 
(von-Neumann) entanglement entropy of a subregion of a free 
quantum field in flat space time, on the other hand for the 
black hole entropy. This intriguing connection was first sug- 
gested and explored in refs.^^ ^''^ and extended in refs."*^ '^^ 
Indeed, there are physical arguments that make the reduc- 
tion of the situation of having a scalar field in a static spheri- 
cally symmetric space-time to a scalar field in flat space time 
plausible^' . The exact status of the relationship between these 
quantities (or to what extent they are related by originating 



V. AREA LAWS FOR CLASSICAL SYSTEMS AND FOR 
TOTAL CORRELATIONS 

A. Classical harmonic systems 

Throughout this article, we have been concerned with quan- 
tum systems on a lattice. What if we have classical systems 
on a lattice, could one still expect an area law to hold? Ob- 
viously, the concept of entanglement is no longer meaningful. 
Also, the Shannon entropy of, say, the marginal distribution of 
a distinguished region / would not quantify correlations in a 
reasonable fashion. What is worse, in case of harmonic classi- 
cal systems on a lattice, when thinking in terms of phase space 
cells, this quantity is burdened with the usual Gibbs para- 
dox. However, it does make perfect sense to talk about classi- 
cal correlations in classical systems, the appropriate quantity 
grasping such correlations being the mutual information: 

Given a probability distribution p on the lattice L, one can 
quantify the correlations between the marginals with respect 
to a distinguished region / and its complement O by means 
of the mutual information. It tells us how much information 
can be obtained on O from measurements in /, and equally 
on / by measurements in O. This quantity enjoyes a number 
of very natural properties. The mutual information is always 
positive — there can be no negative correlations — and will van- 
ish exactly if the probability distribution factorizes, in which 
case one can not learn anything about O from /. Given the 
marginals pj and po of the probability distribution p on / and 
O, respectively, the mutual information is defined as 

I{I:0) = S{pi) + S{po)-S{p), (35) 



For a short general review on this connection see e.g., ref." , for a calcula- 
tion of the one-loop correction to the Bekenstein-Hawking entropy in the 
presence of matter fields and its relationship to the geometric entropy see 
j.gj^j2lo] for an entang lement-based view of the Bekenstein-Hawking 

entropy law see refspma. 
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where here S{p) = — J^j Pj log2 (Pj ) '^he standard informa- 
tion theoretical Shannon entropy. It is noteworthy that the mu- 
tual information does not suffer from the Gibbs paradoxon as 
will be shown below. How does the mutual information scale 
with the size of a region in case of a harmonic coupled classi- 
cal system? The subsequent statement clarifies this situation: 
Consider a classical harmonic lattice system, with Hamilto- 



H = 



■j,keL 



(36) 



where now x = (xi, . . . , XAf) and p = (pi, . . . ,pAr) are 
the vectors of classical position and momentum variables 
of classical oscillators arranged on a cubic lattice L = 
{1, . . . ,iV}^^. The phase space coordinates are then ^ = 
{x,p). The matrix V £ rI-^I^^I^I with a finite-ranged interac- 
tion defines the interaction. 

The state of the system is defined by the phase space den- 
sity, so a classical distribution p : R+. For any 
non-zero inverse temperature /3 > 0, this phase space distri- 
bution is nothing but 



d£,e 



To define the mutual information, following the standard pro- 
cedure, we split the phase space into cubic cells each with a 
volume h?^ , with h > Q being some constant. From the 
phase space space density, we can then identify a discrete 
probability distribution, from an average of the phase space 
density over these cells, pj = /^^^ d^p{£) for j e L. The dis- 
crete classical entropy is then defined as the Shannon entropy 
of this probability distribution as 



Sc{h) = -^Pjlog2(pj). 



We now come back to the situation of having a lattice sys- 
tem with an interior / and an exterior O. The respective 
discrete classical entropies are defined as Si{h) and So{h). 
Obviously, the values of these entropies will depend on the 
choice of h, and in the limit — > 0, they will diverge, log- 
arithmically in h. This is a familiar observation in classi- 
cal statistical physics, the divergence being resolved in the 
third law of thermodynamics. Here, we are, however, inter- 
ested in classical correlations, as being quantified in terms 
of the mutual information which in the limit of h ^ Q 
is well-defined. Hence we can define the classical mutual 
information of a harmonic lattice system as /(/ : O) = 
lim^_).o {Si{h) + So(h) — Sc{h)). We are now in the posi- 
tion to state the area theorem for classical harmonic systemJ^ 

Theorem 16 (Correlations in classical harmonic systems) 

Consider a harmonic lattice system with Hamiltonian as in 
eq. (36\ on a general lattice G — (L, E). Then the classical 
mutual information I{I : O) of the Gibbs state at some 
inverse temperature /3 > satisfies an area law, 

I{I : O) = 0{s{I)). 



The interesting aspect of this prooF^ is that it relates this 
question of the classical mutual information to a quantity that 
arises in the quantum case in case where the coupling matrix 
Vx is replaced by V^, and is hence a simple corollary of ear- 
lier results on quantum systems, now with a coupling that is 
replaced by the squared coupling matrix. Hence, a ''quantum 
proof can be applied to establish a statement on classical lat- 
tice systems. The lesson to learn is that whenever one has 
local interactions — even in classical systems — one should not 
be too surprised if this manifests itself in an area law in the 
correlations. 



B. Classical correlations quantum spin models 

The situation is even simpler for finite-dimensional con- 
stituents. Indeed, quite in contrast to the overburdening tech- 
nicalities that render the question of area laws in higher- 
dimensional quantum systems at zero temperature so difficult, 
the situation can here be clarified with hardly any mathemat- 
ics at all: An elegant, but simple argument shows that total 
correlations in quantum (and classical) systems at non-zero 
temperatures always satisfy an area law. This is a statement 
on correlations — not en tanglement, in contrast to the discus- 
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some /3 > for classical or quantum systems^. The relevant 
quantity grasping correlations is again the mutual information 



I{I:0)^S{pi) + S{po)~S{p), 



(37) 



where S stands either for the von-Neumann quantum entropy, 
or for the Shannon entropy of the probability distribution. The 
classical variant was first discussed in ref.^^, the quantum ver- 
sion in refs.^^-i^. Ref."*^ introduces this quantity to avoid di- 
vergencies of the entanglement entropy in quantum field the- 
ory: In a similar fashion as above, regulators will in fact can- 
cel each other, and the familiar ultraviolet divergence in the 
quantum field limit disappears. 

Interestingly, a general statement on the scaling of corre- 
lations at non-zero temperature in terms of eq. ( [37] i can be 
derived which holds for any spin model with local dimension 
d (see page 295 of ref.'"' and the subsequent ref.lSll): 

Theorem 17 (Classical correlations at non-zero temperature) 

Consider a classical or a quantum system with finite local 
dimension d defined on a translation-invariant lattice 
G — (L, E). Consider the Gibbs state at some inverse 
temperature /3 > of a local Hamiltonian H with two-site 
interactions. In the classical case, where each of the lattice 
sites corresponds to a spin with configuration space TL^, 



l{l:0)<W)\\og{d). 



(38) 



For a quantum system with local Hilbert spaces (C^, the mu- 
tual information satisfies the area law 



I{I:0)<m\W)l 



(39) 



where \\h\\ is the largest eigenvalue of all Hamiltonians across 
the boundary of I and O. 
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This statement is valid in remarkable generality, given the 
simplicity of the argument. We will focus on quantum systems 
in the following. One can write the Hamiltonian H having 
two-site interactions as H = Hj + Hg + Hq, where Hj and 
Ho collect all interaction terms within the regions, whereas 
Hg stands for terms connecting the two regions. The Gibbs 
state pp for some inverse temperature (3 > minimizes the 
free energy F{p) = tr[_ffp] — S{p)/f3. Clearly, therefore, 

Fipp)<F{pi(dpo), 

from which /(/ : O) < l3vi:[Hg{pi®po — pp)] is obtained. As 
the right hand side depends only on terms coupling the inside 
to the outside, i.e surface terms, eq. ( [39] ) follows straightfor- 
wardly. A naive limit /3 — > oo will not yield an area law for 



zero temperature, as the right hand side of eq. ( 39 1 then clearly 



diverges, but for any finite temperature, one obtains a bound. 



VI. CONNECTION TO SIMULATABILITY 

There is an intimate connection between area laws for the 
entanglement entropy and questions of the simulatability of 
quantum many-body systems. The fact that there is "little en- 
tanglement" in a system that satisfies an area law is at the core 
of the functioning of so powerful numerical techniques as the 
density-matrix renormalization group (DMRG) methods. To 
describe the large research field of numerical simulation us- 
ing DMRG-type methods would be beyond the scope of the 
present review. Instead, we will concentrate on the direct rela- 
tionship between the "effective degrees of freedom" that must 
be considered when classically describing quantum systems. 



A. Numerical simulations with the density-matrix 
renormalization group method 

This connection is particularly clear in one-dimensional 
systems, that is for quantum spin chains. Indeed, one can say 
that the fact that ground states of gapped systems satisfy an 
area law — and to a lesser extent that critical systems merely 
have a logarithmic divergence of the entanglement entropy — 
is responsible for the success of the density-matrix renormal- 
ization approach. Matrix-product states also satisfy a one- 
dimensional area law. As MPS are underlying the DMRG 
approach this suggest that the entanglement content of a state 
and the best possible performance of a DMRG approach can 
be intimately linked. 

Historically, DMRG was born out of an idea of renormal- 
ization, where one iteratively identifies the relevant degrees 
of freedom, grasping the essential physics of the problem, 
when going from one step of the procedure to the next one. 
This general idea goes back to the real-space renormaliza- 
tion group approach, presented in ref .1^2^ in the mid 1970ies. 
This approach was particularly successful in the numerical as- 
sessment of the Kondo problem, whereas for other problems, 
results were not quite what was hoped for. The birth of the 
DMRG approach as such was related to a clear analysis of the 
strengths and weaknesses of the real-space renormalization 



group approach to study the low-energy properties of quantum 
many-body systemJSl! Ref.l222| seen as the manuscript in 
which the DMRG method has actually been introduced. Since 
then, this method has seen a standard method in the numeri- 
cal study of strongly correlated quantum many-body systems. 
For a recent review, see ref.CSS]. 

Initially, the formulation of DMRG was based on the above 
renormalization idea. However, in the following years it be- 
came clear that DMRG generates matrix-product states, an in- 
sight that has been reported in rsiW for the thermodynami- 
cal limit of DMRG, and in ref.l^for finite-size DMRG meth- 
ods with the latter placing a particular emphasis on exploit- 
ing a rotational symmetry in variational approaches. Ref.LHSI 
gives a relatively early exhaustive overview over variational 
ansatzes with matrix-product states and the relationship with 
the DMRG idea. Ref.'^^ already hinted at the possibility for 
treating period boundary conditions in the MPS picture but 
chose translation invariant matrices. Ref.^'^ relaxed this con- 
straint to demonstrate that a suitable formulation significantly 
outperforms standard DMRG for periodic boundary condi- 
tions in terms of memory requirements. 

Hence, DMRG — in its several variants — can be seen as a 
variational method, where the optimization problem 



minimize (V'|i?|V'')j 
subject to e (C*)*^^, 



(40) 



impractical already because of its exponentially large feasible 
set, is replaced by a variant of an optimization problem over a 
polynomially large set 

minimize (41) 
subject to 1?/;) G (C'')'^^ is an MPS vector of dimension D. 

In this variant, or — more accurately — in each of these variants 
one does not attempt in one go to identify the global optimum, 
but rather effectively iteratively solves for the local matrices 
involved. Such an iteration will then certainly converge (albeit 
strictly speaking not necessarily to the global minimum).'^ 



B. Approximation of states with matrix-product states 

Any such method, can then only be as good as the best 
possible MPS can approximate the true ground state at hand. 
This, in fact, is related to the entanglement content, in that it 
matters whether or not the true ground state satisfies an area 
law or not. In the light of previous discussions, this connec- 
tion is not that surprising any more: After all MPS satisfy an 
area law for the entanglement entropy. Hence, one aims at ap- 
proximating ground states with states that have in this sense 
little entanglement, and those states can be well approximated 
by MPS that satisfy an area law in the first place. 

This connection has been hinted at already in the first work 
on DMRCpa where the spectrum of the half chain has been 



For mixed state simulations, see refslMni. 
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considered and put into relationship with the "truncation er- 
ror" in DMRG. This is a key figure of merit of the quality of 
an approximation in a step, so unity minus the weight of those 
terms being kept in a step of the iteration. 

This connection between the decay of spectral values of 
half chains, the more rapid the decay the better can DMRG 
perform, has been made more precise and fleshed out in 
j-gflI22l In ref.G^ the relationship to criticality in this con- 
text has been emphasized. Ref.^^^ is a short review on this 
question. In more recent quantitative approaches, the optimal 
approximation that can possibly be obtained by a MPS of a 
given D is considered. Let us denote with Hn ~ ((C^)®^ 
the Hilbert space of a quantum chain of length A^. MPS are 
considered as defined in eq. ( [29] l for open boundary condi- 
tions. Given a family {\iPn)}n of state vectors, it is said 
that it can be approximated efficiently by MPS if for ev- 
ery 6 > there exists a sequence \ipN,D{N)) of MPS with 
D{N) = 0{po\yg{N)) such that 

II |V'Jv)('0Af| - \lpN,D{N)){'lpN,DiN)\ 111 < S, 

where || .|| i denotes the usual trace-norm. In contrast, it is said 
that this sequence cannot be approximated efficiently by MPS 
if there exists some S > such that no sequence of MPS with 
D{N) growing polynomially can approximate up to a 

small error S in trace -norm^"-': 

Theorem 18 (Approximatability with MPS) Consider 
sequences of state vectors {\'ipN)}N G 'Hn of a quantum 
chain of length N, and denote as before the reduced state of 
a block / = {1, . . . , n} q/ length n with pj. If the sequence 
of pj satisfies an area law for a Renyi entropy S^for a < 1, 

SM = 0(1), 

then the sequence {\4'n)}n is efficiently approximable by 
MPS. In contrast, if the von-Neumann entropy Si{pj) = 
r2(n), so grows at least linearly with the block size, then it 
cannot be approximated efficiently by MPS. This means that 
states satisfying a volume law cannot be approximated. The 
same holds true if any Renyi entropy S^for some a > 1 grows 
at least as Sa{pi) — ^{n'^) for some k < 1. Otherwise, the 
connection is undetermined, in that examples both for approx- 
imable and inapproximable states can be found. 

This statement clarifies the connection between the entan- 
glement content and the possibility of describing states with 
matrix-product states. The validity of an area law implies that 
there is sufficiently little entanglement in the state such that 
an economical description in terms of matrix-product states is 
possible. The enormous success of DMRG is related to the 
fact that gapped systems satisfy an area law. Even if the sys- 
tem is critical, the logarithmic divergence still allows for a 
relatively economical description in terms of matrix-product 
states. The fact that Renyi-entropies for a smaller than or 
larger than unity feature here may be seen rather as a technical 
detail. The general message is clear: The area-like entangle- 
ment scaling, with or without small corrections, allows for an 
efficient approximation in D for matrix-product states. 



To reiterate the point made in Subsection IV.G Quenched, 
non-equilibrium systems can indeed fall exactly into the cate- 
gory of having an effectively linearly growing block entropy, 
so are characterized by a volume law for the entanglement en- 
tropy. More precisely, we face the interesting situation that 
while for each time, we have an area law in n, the constant in 
the upper bound grows in time such that for a suitable choice 
for the sub-block, one arrives effectively at a volume law, as 



made precise in Theorem 15 This has severe pract ical impli- 
cations: For small times, t-DMRG^-^^^^^ , the time- 
dependent version of DMRG, can very accurately keep track 
of the dynamics of the system. This is a variant in which one 
essentially makes a Lie-Trotter approximation of the time evo- 
lution operator, and then approximates in each time step the 
resulting state vector by an MPS, going back to ref.^^'. The 
functioning of this algorithm can essentially be traced back to 
the observation that an arbitrarily good approximation to the 
propagator can be established with polynomial computational 
resources in the system sizd^^. In time, however, one will 
eventually encounter typically an exponential increase in the 
number of degrees of freedom to be kept in order to faithfully 
describe the state. This eventually limits the time up to which 
one can numerically simulate time evolution using a variant of 
DMRG. The increase in the entanglement content also even- 
tually limits classical simulations of quantum adiabatic al- 
gorithms based on MPS, which nevertheless perform often 
impressively well (for a careful numerical analysis, see, e.g., 
ref ^). It is interesting to note, however, that this complexity 
does not necessarily translate in the difficulty of following the 
time-evolution of specific observables when evolving them in 
the Heisenberg picture using t-DMRG. Then, in some cases 
the Heisenberg time evolution can be carried out exact ly for 
finite bond dimension and arbitrary long time^^lIIlIEIIl 

There are numerical simulation methods that allow for the 
simulation of certain quantum states that do not satisfy an area 
law. MERA already allows for a logarithmic divergence of the 
entanglement entropy in one-dimensional systems. Weighted 
graphs state based approached and its 1-D variant, the renor- 
malization algorithm with graph enhancement^ can cope 
with instances of volume laws for the entanglement entropy, 
the latter in 1-D the former in arbitrary spatial dimensions. 
Early work on the simulation of a particular kind of discrete 
time evolution, namely the application of random unitary cir- 
cuits, suggests that this may be a promising approach for the 
"efficient simulation of quantum many-body systems beyond 
area laws". 

We end this subsection with a note rather from the com- 
puter science than from the physics perspective: The fact that 
a true ground state is well-approximated by an MPS does, 
strictly speaking, not necessarily mean that DMRG will also 
efficiently find this best approximation. In practice, DMRG 
works well, and it typically produces good and reasonable re- 
sults. It is remarkable how well this approximation is found 
in the iterative scheme as being pursued by any DMRG al- 
gorithm: After all, the full problem eq. (41 1 is a non-convex 



polynomial global optimization problem of very high order 
((V'lff IV') is of degree N"^ in D). Still, by local variations 
and sweeping one achieves very good results. The ultimate 
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reason for this impressive performance is yet to be ultimately 
understood. 

Having said that, the worst case complexity of the problem 
of finding the best approximation can be computationally dif- 
ficult in the sense of computer science. In fact, the class of 
problem of keeping some matrices fixed and varying over a 
finite subset has in worst case instances that are NP-hara^^. 
In non-translation invariant settings, one even finds that if one 
could efficiently identify the best possible MPS approxima- 
tion, one could solve efficiently NP-hard problems^**. Even 
more strongly put, the problem of approximating the ground 
state energy of a system composed of a chain of quantum sys- 
tems is QMA-complet^. 

This should be seen as a warning sign: The functioning of 
variational algorithms such as DMRG is essentially based on 
heuristics, and in worst case one can encounter hard problems. 
The energy landscape is then so rugged that one gets stuck in 
local optima. Still, while it is important to acknowledge that 
DMRG is strictly speaking not certifiable, it is still true that 
it works very well in practice and is one of the pillars of the 
numerical assessment of strongly correlated systems in 1-D. 



C. Implications on higher-dimensional simulations 

For higher dimensional systems, tensor-product states or 
PEPS, as well as those of MERA, satisfy area laws, as has 



been discussed in Subsection IV.F This fact suggests that 
when minimizing {il]\H\ip) for an x A^-lattice subject to 
IV') e (C'')®^' being a PEPS or MERA described by polyno- 
mially many real parameters, one encounters a good approxi- 
mation whenever the system at hand already satisfies an area 
law. In the light of the fact that even critical two-dimensional 
systems can satisfy an area law, this would mean that they can 
be well-described by PEPS or MERA described by relatively 
few parameters. Numeri cal work in c ase of PEPS indicates 
that this is indeed the casJSIEllIIIIIIIl gyen in the t hermody- 
namical limit of an infinite systerrF^or for fermionJ '^l'"^ l '^° [ 

A rigorous result similar to TheoremfTS] yet, is still lacking 
for PEPS or MERA. The intuition developed so far, however, 
is in one way or the other quite certainly right: Whenever an 
area law is satisfied, PEPS with small bond dimension should 
give rise to a reasonably good approximation. Here, subtle 
aspects are rather connected to the fact that the exact contrac- 
tion of the tensor networks of PEPS, and hence the compu- 
tation of expectation values, is inefficient, and that approx- 
imate contractions have to be employed. Suitable subsets, 
such as the class of string states, can always be efficiently 
contracted, giving rise to v ery promising variational sets in 
higher-dimensional systemJ^H The method in ref.SS also 
gives rise to certifiable approximations of 2-D ground states 
for a class of models, exploiting quasi-adiabatic evolutions. 

As before, one has to distinguish the variational set as 
such from the computational method of varying over this set. 
Usually, one has to find practical and heuristically suitable 
methods of solving a global optimization problem over many 
variables. Several strategies may be followed when varying 
over suitable sets to simulate higher-dimensional strongly cor- 



related systems: One may use local variations such as in 
DMRG, imaginary time evolution, or flow method^, mak- 
ing use of gradient flow and optimal control ideas to vary over 
the manifold of unitary gates that describe the variational set 
of states at hand. For MERA, the same intuition should hold 
true. Here, first approache s implem ented have been focused 
on one-dimensional system^^-^'™''^^', but the ideas are also ap- 
plicable in higher dimensions'"* ^- 



VII. PERSPECTIVES 

In this Colloquium, we have presented the state of affairs in 
the study of area laws for entanglement entropies. As has been 
pointed out above, this research field is presently enjoying a 
lot of attention, for a number of reasons and motivations. Yet, 
needless to say, there are numerous open questions that are 
to be studied, of which we mention a few to highlight further 
perspectives: 

• Can one prove that gapped higher-dimensional general 
local lattice models always satisfy an area law? 

• In higher dimensional systems, critical systems can 
both satisfy and violate an area law. What are further 
conditions to ensure that critical systems satisfy an area 
law? What is the exact role of the Fermi surface in the 
study of area laws in fermionic critical models? 

• Can one compute scaling laws for the mutual informa- 
tion for quasi-free systems? 

• For what 1-D models beyond quasi-free and conformal 
settings can one find rigorous expressions for the entan- 
glement entropy? 

• Under what precise conditions do quenched disordered 
local models lead to having "less entanglement"? 

• What are the further perspectives of using conformal 
methods for systems with more than one spatial dimen- 
sion? 

• Can the link between the Bekenstein formula in the AdS 
context and the scaling of geometric entropies in con- 
formal field theories be sharpened? 

• To what extent is having a positive topological entropy 
and encountering topological order one to one? 

• How can the relationship between satisfying an area law 
and the efficient approximation of ground states with 
PEPS be rigorously established? 

• What efficiently describable states satisfy an area law, 
such that one can still efficiently compute local proper- 
ties? 

• Are there further instances for 1-D systems satisfying 
an area law that allow for certifiable approximations in 
terms of matrix-product states? 
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These questions only touch upon the various perspectives that 
open up in this context. The quantitative study of a research 
area that could be called "Hamiltonian complexity"^" is just 
beginning to emerge. The puzzle of how complex quantum 
many-body systems are, and how many effective degrees of 
freedom are exploited by nature, is still one of the intriguing 
topics in the study of interacting quantum systems. 
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IX. APPENDIX: FISHER-HARTWIG THEOREM 

In this appendix we briefly present an important techni- 
cal result concerning the asymptotic behavior of Toeplitz 
matriced^. 

Lemma 1 (Fisher-Hartwig) Consider a sequence of n x n 
Toeplitz matrices {T„}„ with entries {Tn)ij — (T„)i_j, 

generated by g : [0, 27r) — >■ (D. Let g be of the form 

R 

r=l 

with tp{(j)) = e-^0(^-'l>)^ Ua = {2- 2cos((/)))°, Re{a) > 
— 1/2, and b : [0, 27r) — > (D a s mooth no n-vanishing function 
with winding number zero. 

and \Reil3r)\ < 1/2 or R = 1, a = \Re{P)\ < 5/2, the 
asymptotic behavior of the determinant ofTn is given by 

y det{Tn) _ 



This terni has been coined by B.M. Terhal. 



where E = 0(1) in n and 

G = exp(^i-^ d01n(6(0)) 
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